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_…— Deze hoofdstukken, grootendeels herziene overdrukken van sommige 
} mijner opstellen in buiten- en binnenlandsche tijdschriften (l’Eseignement 
_J. Mathématique; Jahresbericht der Deutschen Mathematiker-Vereinigung; 
… Astronomische Nachrichten; Wiskundig Tijdschrift; Nieuw Tijdschrift voor 
4 Wiskunde), zijn bestemd voor hen, die Wiskundige acten trachten te 
behalen, — met name Ki. 

Bovendien voor de studenten ín exacte Wetenschappen aan onze 

59 Universiteiten en aan de Delftsche Hoogeschool. 

| Ook aan onze Militaire Academie en aan de Landbouwschool te 
Wageningen, waar de Wiskunde naar verdienste wordt behartigd, kunnen 
zij wellicht menigeen van dienst zijn. 

In het algemeen zullen zij, naar ik hoop, ieder, die de Mathesis met 
ambitie beoefent, belang kunnen inboezemen. 

De titel dezer Ge opstellen moge hier even nader worden 
toegelicht. 

„Wegens den reuzenarbeid, Ean groote genieën reeds gedurende eenige 
eeuwen in de Hoogere Algebra verricht, is het practisch ongeveer on- 
mogelijk nog zuiver maagdelijk gebied te vinden of geheel nieuwe wegen 
aan te leggen, die nergens reeds gebaande paden zouden kruisen of volgen. 

Onder „Nieuwe Methoden’ wordt daarom de beoordeelaar verzocht te 
willen verstaan: nieuwe aanwendingen, wijzigingen en verbeteringen en 
ook nieuwe combinaties van methoden op het gebied der hoogere- 
machts-vergelijkingen. 

Mijn werk is echter, — hetgeen naar ik hoop bij nader onderzoek zal 
blijken, — oorspronkelijk, — voor zoover dit eenigszins mogelijk is. 

Hadden voor mij andere wiskundigen reeds eenig denkbeeld over 
sommige quaesties of methoden geopperd, dan is daarvan rekenschap 
afgelegd. 

Ik hoop werkelijk hier en daar eenig Hate nieuws te hebben gebracht. 

De inhoud van dit werkje spreke verder voor zich zelf. Alleen meen 

tr ik hier, — ter voorkoming van een verkeerden indruk, — wel te moeten 
„opmerken, dat de door mij gegeven ongewone, geheel elementaire atlei- 
“ding eener implicite ontwikkeling niet aldus is behandeld wegens onbe- 
Skendheid met desbetreffende hoogere theorieën, maar wel om deze 
_rquaestie en hare zeer nuttige toepassingen binnen het terrein der studie 
evoor Kr te brengen. 
3 De met zorg geconstrueerde grafieken, die als Aanhangsel zijn toege- 
‚a voegd, zijn behalve ter illustratie van de beschouwingen over de 
“Kepler'sche Vergelijking ook bestemd voor het werkelijk in praktijk 
brengen der voorgestelde graphische benaderingsmêthoden. 

Aan alle astronomische observatoria kunnen zij, naar ik vast vertrouw, 


goeden dienst bewijzen. 
Dr. B. GONGGRIJP. 
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A. Grafiek van de Sinusoide ter (benaderde) graphische oplossing 
der Kepler'sche Vergelijking. 

B. Gewijzigde Sinusoide. 

C. Drie systemen rechte lijnen. 


L. OPLOSSING VAN HOOGERE-MACHTS-VERGELIJKINGEN 
DOOR MIDDEL VAN KETTINGBREUKEN. 


Het principe dezer methode is reeds van Lagrange afkomstig 
en uiterst eenvoudig. 
Zij x een reëele wortel van de vergelijking: 
Ln 0 dn Antidns Aril 
Wij denken het aantal geheelen, a, van dezen wortel bepaald 


en stellen: 
Nm (l + Eee 
Xi 
In de nieuwe vergelijking, waarin de onbekende nu x, is, 
stellen wij verder: | 
l 





ar One 
waaruit volgt: 
l 
ond erna 
a, en PR 
X3 


bed 


a, is nu het aantal geheelen van x, in de nieuwe vergelijking. 

Op deze wijze kan men voortgaan. Het ís duidelijk, dat zoo- 
doende x in een kettingbreuk ontwikkeld wordt. 

Dit beginsel is zoo eenvoudig en de ontwikkeling van een te 
benaderen getal in een kettingbreuk heeft zulke voordeelen }), dat 
men zich met verbazing afvraagt hoe het komt, dat de methode 
geheel veronachtzaamd is en in de leerboeken over Hoogere 
Algebra alleen maar genoemd wordt, gewoonlijk met de opmerking, 
dat de z.g. „Horner'sche methode” de voorkeur verdient. 

De verklaring van deze verwaarloozing moet geloof ik gezocht 
worden in het feit, dat Lagrange er niet in geslaagd is de noodige 
vereenvoudigingen aan de voor zijn methode vereischte bewer- 
kingen aan te brengen. Daarentegen werd door invoering van 
den Horner'schen Algorithmus bij de Newton'sche methode deze 
laatste zoozeer verbeterd, dat zij de andere geheel in de schaduw 


| É 
1) Men denke hierbij aan de zeer bijzondere eigenschap: Als — een „naderende breuk” 
n 
is van een kettingbreuk-ontwikkeling, waardoor een getal G wordt voorgesteld, dan 
3 1 
verschilt 3 minder dan —- van G. Tengevolge van deze eigenschap behoeft zulk een 
n° 


naderende breuk slechts een vrij bescheiden noemer te hebben om het getal G met 
groote nauwkeurigheid te kunnen voorstellen. Een beroemd voorbeeld is de benadering 
‘van Metius voor het getal zr, nl.  — 355/,, 3. 


} 


stelde. Zwaar zal ook de omstandigheid gewogen hebben, dat door 
de oplossingsmethode van Newton spoedig elk volgend cijfer van 
den wortel zonder eenigen twijfel kon worden gevonden, terwijl 
bij die van Lagrange het bepalen van elk volgend wijzergetal van 
de kettingbreuk een afzonderlijk onderzoek schijnt te vereischen. 

In de volgende bladzijden hoop ik aan te toonen, dat de methode 
van Lagrange zeer belangrijk vereenvoudigd kan worden o.a. door 
ook deze te combineeren met den Algorithmus van Horner; dat 
zij dan gerust naast die van Newton mag gesteld worden en 
dikwijls zelfs boven deze de voorkeur verdient. 

Allereerst dient te worden aangetoond, dat men steeds op 
eenvoudige wijze de grootheden a, A, dy enz. kan vinden, die 
bij een zelfden wortel behooren. 

Denken wij voorloopig alleen reëele wortels en de methode 
toegepast op den grootsten positieven wortel. 


Wij stellen: geese of 
Xi Xi 


Wij verminderen dus alle wortels eerst met a. De grootste wortel 
der oorspronkelijke vergelijking gaat dan over in een positieven 
wortel der nieuwe verg. met getallenwaarde kleiner dan 1; de 


l 
overige wortels worden negatief. Vervangt men nu x door Se 
1 


m.a.w. keert men nu de wortels om, dan ontstaat een verg., die 
slechts één positieven wortel zal hebben; zooals van zelf spreekt 
zal deze wortel dan grooter dan 1 zijn. Men kan dan onmiddellijk 
dezelfde bewerking herhalen. 

De methode geldt echter voor elken positieven en ook voor 
elken negatieven wortel mits men bij enkele van de eerste trans- 
formaties eenige omzichtigheid aanwendt. 

Nemen wij een concreet geval. Gesteld, een vergelijking heeft 
3 positieve wortels, waarvan wij voorloopig aannemen, dat allen 
een verschillend aantal geheelen bezitten. 

Zij nu a het aantal geheelen van den kleinsten wortel. Stellen 
wij eerst: x=aty, of £—a =y, dan zal de resulteerende 
vergelijking in y nog drie positieve wortels hebben, maar een 
er van ligt tusschen O en l en de andere twee zijn grooter dan 1. 

Keert men dan de worters om door de substitutie y = Ee dan 
zal de vergelijking in x, nog slechts 1 wortel > 1 kunnen hebben; 


d 


de beide andere zijn dan < 1. Men ís nu teruggekomen op het 
eerste geval. 
Wanneer er twee wortels tusschen dezelfde twee geheele ge- 
tallen gelegen zijn, zal men ze door de substitutie 
Sleen Ee 
X% | 
en misschien zelfs door enkele volgende niet terstond kunnen 


scheiden, maar weldra zal deze scheiding toch moeten plaats 
vinden, — tenzij de twee gedachte wortels gelijk zijn. Ook in 
dit geval kan, zooals wij later zullen toelichten, de methode 
blijven dienst doen op geheel analoge wijze als die van Newton. 
De bewerkingen, bij de boven aangegeven methode vereischt, 
schijnen op het eerste gezicht vrij omslachtig te zijn en door het 
telkens bepalen van het aantal geheelen van een wortel eenigermate 
op den tast verricht te moeten worden. Wij zullen nu echter de 
vereenvoudigingen bespreken. 
Wij stellen eerst „=aty of x—a=y 
en bepalen door den Algorithmus van Horner de coëfficienten 
der vergelijking in Xx — 4. 


. k l 
Hierna zouden wij y — Pe moeten stellen en de nieuwe verge- 
1 


lijking, gerangschikt naar x,, moeten opschrijven. 

Hiertoe is echter niets anders noodig dan de volgorde der 
coëfficienten der vergelijking in y om te keeren, m.a.w. ze van 
rechts naar links te lezen. 

Daarna bepalen wij het aantal geheelen van x; zij dit a, 
(Over vereenvoudigingen, die hierbij optreden, spreken wij nd 

Nu moet allereerst weer gesteld worden: 

| Kie MOL Wid 

De wortels der vorige vergelijking, (waarvan de coëfficienten 
van rechts naar links gelezen worden), moeten nu met a, ver- 
minderd worden. 

Men behoeft nu de coëfficienten daarvan in het geheel niet 
opnieuw op te schrijven; men kan eenvoudig den Horner’ schen 
Algorithmus nu van rechts naar links toepassen. 

Bij de achtereenvolgende stappen in de benadering past men 
dien Algorithmus nu steeds afwisselend van links naar rechts en 
van rechts naar links toe. 

Beter dan lange uitweidingen zal een voorbeeld dit verduidelijken. 


1 


Voorbeeld. xs HX? — ZX — 1 =0. 

Deze vergelijking heeft slechts één positieven wortel, die bij 
onderzoek tusschen l en 2 blijkt te liggen. Wij hebben dus op 
de gewone wijze de wortels der vergelijking met 1 te verminderen: 


1 1 ES bes dl 
(A) 1 9 0 
ND WOT add: 
1 3 È 
8 ES AEG 
1 
4 


Nu zouden wij de wortels hebben om te keeren, zoodat wij 

verkregen: | 
— 8 —_ 4 — Ì 

(De teekens zijn nu ook omgekeerd). 

Wij mogen nu van de wortels der vergelijking, welker coëffi- 
cienten door deze rij worden voorgesteld, zooveel aftrekken, dat 
de laatste coëfficienten nog juist zijn teeken behoudt, maar het 
tegenstelde teeken zou verkrijgen, als men een eenheid meer van 
de wortels ging afnemen }) 

Men ziet zeer gemakkelijk, dat 4 het grootste getal is, dat 
onder deze voorwaarden van de wortels kan worden afgetrokken. 

We zullen dus de bewerking hebben: 


eet 
4 4 0 
rent 
4 20 
WADE G20 
4 
RE 


Van de opnieuw te onderzoeken vergelijking zullen de coëffi- 
cienten dan zijn: 
l — 20 —g9 — | 
Nu is het onmiddellijk duidelijk, dat het voor het vinden van 
dit resultaat niet noodig geweest was de teekens van de rij 
— | 5) 4 l 
om te keeren. Men kan evengoed de wortels der vergelijking 


1) Algemeener: zooveel aftrekken, dat de eenige variatie nog juist blijft bestaan, 
maar zou verdwijnen, als met een eenheid meer werd verminderd. 
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met negatieven eersten coêfficient met 4 verminderen en dus de 
volgende bewerking verrichten: 


— | 5) 4 1 
— 4 — 4 0 
EE EDER, POR GOTE 
— d — 20 
En E20 
—4 
0 
Men ziet nu juist de coëfficienten-rij 
| — 20 — 9 — Ì 


te voorschijn komen. 

Welnu, daar dit zoo is, kunnen wij de tweede bewerking 
onmiddellijk aan de bewerking (À) laten aansluiten door den 
Algorithmus van Horner van rechts naar links toe te passen. 

De twee bewerkingen geven dan te zamen het volgende tableau: 


1 1 MN) EE) 
0 1 2 0 
fl 0D 0 EN 

1 3 

So tds 

1 add 

BZ Bart 

de Ee 

BD Rn 

0 Hee 

200 ENE 


Op deze wijze gaat men steeds voort. De Algorithmus wordt 
steeds afwisselend van links naar rechts en dan weer van rechts 
naar links toegepast. | 

Een belangrijke questie is nu nog de bepaling van het aantal 
geheelen, waarmee telkens de wortels verminderd moeten worden. 
Wij zullen aantoonen, dat men zich hieromtrent zeer gemakkelijk 
alle zekerheid kan verschaffen. | 

In ons voorbeeld was er slechts één positieve wortel, nl. tusschen 
l en 2. Wij schreven daarna de vergelijking naar (x — 1) op. 
Deze had dus één positieven wortel tusschen O en 1 en twee 
negatieve wortels, wier volstrekte waarde zeker > 1 was. Door 
ze om te keeren onstond de vergelijking met de coëfficienten-rij: 

— | 5) 4 ie 
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die dus één positieven wortel had > 1 en 2 negatieve, welker 
volstrekte waarde noodzakelijk < 1 is. 

Van deze vergelijking werden de wortels weer verminderd met 
het aantal geheelen van den positieven wortel, zoodat de resul- 
teerende vergelijking weer twee negatieve wortels verkrijgt met 
volstrekte waarde > 1!) en een positieven tusschen OQ en 1. 

En ook deze vergelijking wordt weer omgekeerd, zoodat er 
weer een nieuwe vergelijking ontstaat met één positieven wortel 
> l en twee negatieve, wier getallen-waarde een breuk is tus- 
schen — 1 en 0. 

Dit zal altijd zoo blijven. 

Beschouwen wij nu, om bij ons voorbeeld te blijven, de rij: 

l — 20 — 9 4 
dan weten wij, dat de algebraïsche som der wortels dezer verge- 
lijking gelijk is aan + 20. 

Noemen wij den positieven wortel p en de volstrekte waarden 
der beide negatieve wortels n, en 79, dan moet: 

ph r n3) S= 207 Of: p= 20E AN 

Omdat n‚ en n, beide kleiner dan 1 zijn, zien wij dus, dat: 

pr 20ste Tevens nis EA 

Wij zien, dat het aantal geheelen van p minstens —= 20 moet 
Zijn, maar kleiner dan 22. 

Men zal verder nog het volgende inzien: 

Wanneer het te voren afgetrokken getal (dus het laatst-gevon- 
den wijzergetal der te berekenen kettingbreuk) slechts 2 is, zullen 
NH, en #9 (ontstaan door omkeering van getallen > 2) elk < !/,, 
dus samen < 1 moeten zijn. In dat geval zouden wij hebben: 

JOD 1e 2D eral 
of in het algemeen: 


JAN 
el RET 3 
DIe. OL 


wanneer A, en A, resp. de coëfficienten van de n° en van de 
n— 1° macht der onbekende in de telkens afgeleide vergelijkin- 
gen voorstellen. | 

Als dus het voorafgaande wijzergetal > 2 is, bedraagt het 


aantal geheelen van p hoogstens dat van — ne of l meer; is dat . 
0 


1) Want van een reeds negatief getal, hoewel in volstrekte waarde < 1, wordt 
een getal > 1 afgetrokken. 
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wijzergetal 2 of 1, dan kan daarbij nog één eenheid komen. In 
allen gevalle blijkt echter dit aantal geheelen met genoegzame 
zekerheid bepaald te zijn. In ons geval was het vorige wijzergetal 
reeds 4, zoodat het aantal geheelen van p < 20 + 1/, moest zijn 
en dus onmiddellijk geheel bepaald was. 

Bij vergelijkingen van hoogeren dan den 34e" graad zouden wij 
bijv. moeten stellen: 


pet aan ie or Ha): 


Wanneer dan het voorafgaande. wijzergetal > 4 of —= 4 was, 


zou het aantal geheelen van — 2 hoogstens met de eenheid ver- 
0 


meerderd moeten worden om dat van p te vinden; bij kleinere 
voorafgaande wijzergetallen zou men misschien 2 of 3 eenheden 
meer moeten nemen; in het algemeen zal men echter op het 
eerste gezicht terstond een juiste keus kunnen doen. Bovendien, 
heeft men misschien een eenheid misgetast, dan kan men zonder. 
eenige moeite den Algorithmus nog even toepassen met + 1 of — 1 
al naarmate men het aantal geheelen van p te klein of te groot 
gerekend heeft. Men heeft zulk een tot op 1 eenheid foutieve 
schatting van een cijfer of decimaal ook wel bij de Newton'sche 
methode en zelfs bij de gewone deeling. 

Hieronder hebben wij de nu voldoende duidelijk gemaakte 
methode op het gekozen voorbeeld toegepast. Wij meenen op 
de volgende voordeelen te mogen wijzen: 

1°. Alle bewerkingen geschieden met geheele getallen ; hier- 
door wordt het overbodig om, — zooals bij de methode van 
Horner noodig is, — bij de vermenigvuldigingen en optellingen 
telkens 1, 2, 3 enz. cijfers (van links naar rechts) „uit te schieten”. 

2%. Deze bewerkingen kunnen in het algemeen tot veel kleinere 
getallen (nl. vermenigvuldigtallen) beperkt blijven. 

3%. Wanneer er wel een groot wijzergetal optreedt, waardoor 
bewerkingen met grootere getallen zouden optreden, dan is de 
kettingbreuk ook bijna terstond met voldoende nauwkeurigheid 
bekend. 

40 Reeds van het begin af kan men in het algemeen met 
groote zekerheid het aantal geheelen bepalen, dat van de wortels 
der successievelijk optredende vergelijkingen moet worden afge- 
trokken; dit is bij de methode van Horner niet altijd het geval. 
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5% De deeling, waardoor dit getal gevonden wordt, kan met 
veel eenvoudiger getallen verricht worden. 

6%. Er komt nooit een O in de bewerking, een geval, waarin 
ook bij de methode van Horner afzonderlijk moet worden voor- 
zien en dat daarbij aanleiding kan geven tot vergissingen. 

7%, Teller en noemer der berekende breuk zullen betrekkelijk 
kleine getallen zijn; opdat zij nl. tot op 0,000001 nauwkeurig de 
waarde van den wortel zal voorstellen, behoeft de noemer hoogstens 
slechts uit 4 cijfers te bestaan. 

Men controleere deze opmerkingen aan het geheel uitgewerkte 
voorbeeld hieronder. Als een soort contrôle, dat de bewerkingen 
in voldoend aantal zijn verricht, kan nog het feit dienen, dat na 
het bepalen van een even aantal wijzergetallen de beide uiterste 
en ook de beide binnenste kolommen uit evenveel getallen moeten 
bestaan (dus even lang moeten zijn). 




















l l — 2 —Ì 
0 l 2 0 
l 2 0 — Ì 
196 | 3 — 180 
197 3 3 — 181 
362 Ì —d — 1878 
599 4 —Ì — 2059 
47320 — 4 — 4 — 462 
47879 0 — 5d — 2521 

— 20 — 4 

— 20 —g 

20 0 

0 — 9 Wijzergetallen: 
20 400 1, 4, 20,72, 0e 

20 391 

20 — 362 

40 29 

58 — 362 

98 — 333 

— 666 — 362 

— 068 — 695 

591 69 

23 — 626 

591 1842 

614 1216 

591 — 2059 





Zie volg. blz. 























1205 — 843 
— 843 — 2059 
362 — 2902 
— 9902 — 2059 
— 2540 — 4961 
3354 4884 
814 EEM 
3354 _25008 
4168 24931 
3354 — 95210 
7522 WE 
— 2790 253101 
4732 — 25489 
— 254890 Sn D5 0, 
— 250158 — 50699 


Zooais bekend is, bepaalt men uit boven gevonden wijzerge- 
tallen de naderende breuken van de kettingbreuk aldus: 
En daarna het eerste wijzergetal, 


gedeeld door de eenheid, dus Ee dan vermenigvuldigt men 


Men schrijft voorop altijd 


teller en noemer dezer tweede breuk met het tweede wijzergetal 
en telt bij deze producten den en teller en den vorigen 


noemer op. Dit geeft in ons voorbeeld i ee En Ee 





Vervolgens worden teller en noemer van deze breuk weer ver- 
menigvuldigd met het volgende wijzergetal; beide uitkomsten 
worden weer met den vorigen teller of noemer vermeerderd. Bij 
ons voorbeeld vinden wij dan als volgende naderende breuk 
En = zt Op deze wijze gaat men steeds voort. De 
verdere naderende breuken zullen worden: 

DO A ned 2 0296 63889 
166’ 579’ 745’ 5049 “" 51235: 


De laatste breuk moet reeds tot op minder dan Abed rl of 


25 X 108 
0,00000 00004 nauwkeurig zijn. 

Zonder de geheele bewerking verder voort te zetten ziet men 
gemakkelijk, dat het volgende wijzergetal een 5 zal zijn. De 
825741 
261224 











hierbij behoorende naderende breuk zal zijn: 
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De volstrekte fout in deze laatste breuk moet < Bee, zijn. 

Zij geeft dus zeker in fien decimalen nauwkeurig den positie- 
ven wortel van onze vergelijking. — 

Evenals bij de methode van Horner kan men hier moeite uit- 
sparen door het weglaten van overtollige cijfers. Heeft men nl. 
voor de coëfficiënten der resulteerende vergelijkingen reeds tame- 
lijk groote waarden verkregen, dan kan men alle coëffrcienten 
door 1000, 10000 enz. gedeeld denken en de decimalen wegla- 
ten; deze zullen op de grootte van de nog te volgen wijzerge- 
tallen geen invloed uitoefenen. 

Bij het voorkomen van bijna gelijke wortels, zal elk wijzergetal 


slechts een deel kunnen zijn van wien: Bij twee bijna gelijke 
0 


wortels zal men telkens in de nabijheid van de helft van dat 
bedrag moeten blijven, totdat plotseling de beide wortels der ver- 
kregen vergelijking zich laten scheiden; voor beide kan dan de 
methode afzonderlijk worden voortgezet. 

Dit wordt toegelicht met onderstaand voorbeeld. Te voren zij 
nog opgemerkt, dat het voor ieder, die de Newton'sche methode 
meermalen heeft toegepast, duidelijk zal zijn, dat bij dit voort- 
durend omkeeren der vergelijkingen en het telkens met een aantal 
geheelen verminderen van de wortels tegelijkertijd de bestaanbaar- 
heid der wortels onderzocht wordt en de waarde er van berekend. 

Zij nu op te lossen de vergelijking: 

6x% — 141x + 263 = 0. 

De 2 teekenvariaties houden stand, wanneer men de wortels 

met 2 vermindert. De coëfficienten-rij wordt dan: 
6 36 ee BO 2D. 

Als de wortels met 3 worden verminderd, gaan beide variaties 
echter verloren. Er kunnen dus twee wortels gelegen zijn tusschen 
2 en 3. Deze twee wortels kunnen nu wortels van de nieuwe 
vergelijking worden, maar liggen daar dan tusschen O en 1. Wij 
keeren de wortels om en lezen de coëfficiënten van rechts naar 
links. De bedoelde kleine positieve wortels gaan dan over in 
andere, die weer positief zijn, maar nu > |. 

Wanneer nu ook deze omgekeerde waarden nog in elkaars 


nabijheid liggen, zal men uit p = e <t zeker. een‚te Pan 


pl 


_ waarde voor p vinden, als nu n de volstrekte waarde van den 
eenigen negatieven wortel voorstelt, die onze verg. hebben kan 
en ook zeker heeft. 

Bij twee positieve wortels zal men hebben: 


69 69 
ps, Epa = ig tr en dus p, +p, < 55 + Yn 


want het vorige wijzergetal was een 2. 

Wel volgt hieruit, dat het Kleinste wijzergetal p, zeker 
83,5 | 
2 X 29 

Verminder nu de wortels met 1 (van rechts naar links!), 

De nieuwe rij van coëfficienten wordt: (te lezen van rechts 

naar links!) 


< zal zijn, dus —= l zal moeten zijn. 


ENZ — 15 18 IN, 

Er zijn nog 2 variaties aanwezig; deze zouden echter verloren 
gaan, als men met 2 verminderde (of wanneer men nog eens 
met 1 verminderde). Van beide wortels hebben wij nu de wijzer- 
getallen 2, 1; zij zijn nog niet gescheiden. 

Wij keeren de vergelijking om en lezen dus bovenstaande 
getallen weer van links naar rechts. 

Voor het nieuwe wijzergetal hebben wij: 


lant ed Ee +1, (tvorig wijzergetal was = 1) 


Pr 

AUS pi 4 
Verminder nu de wortels met 4. Er komt: 
2 9 nde sr 


Er is 1 variatie verloren gegaan; de wortels zijn gescheiden 
en hunne bestaanbaarheid is aangetoond }). 

Voor p = 3 zouden er nog 2 variaties geweest zijn. De kleinste 
wortel der laatste: vergelijking heeft dus 3 tot aantal geheelen; 
voor den grootsten hebben wij: 


17 1 
BSD 30 


b Bij toepassing van de Newton-Hornersche methode zou men na bepaling der 
coëfficienten-rij 6 36 — 69 29 
door probeeren moeten vinden, dat bij vermindering der wortels met 0,7 de twee 
variaties nog bleven bestaan, maar tot één gereduceerd werden, als de wortels met 0,8 
werden verminderd. Men vindt dan nl.: 

Be 6 48,6 — 9,78 0,398 en 

20, 6 50,4 + 0,12 — 0,088. 
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Neem dan p = 5 (of verminder nog eens met 1). Er komt de rij: 
2 15 18 — 6. 
Voor p = 6 zou de laatste variatie weer verloren gaan. 
Wij hebben nu voor den grootsten positieven wortel voorloopig 
de wijzergetallen: LEDS 
en als coëfficientenrij om den algorithmus voort te zetten (van 
rechts naar links): 


Ee C:end : 
Voor den Kleinsten positieven wortel vonden wij reeds als 
wijzergetallen: Led Eh 


de coëfficienten, voor het voortzetten van den algorithmus benoo- 
digd, vinden wij bijv. door de wortels van de vergelijking met 
de coëfficienten 2, 9, — 6 en — ll weer te verminderen met 


il; Erkomt dan: 
9 8 18 3 


(Ook te lezen van rechts naar links). 
De verdere bewerkingen verloopen geheel normaal. 
A. Ontwikkeling van den grootsten positieven wortel 


5 15 18 nd: 
45 0 Eng 18 
47 15 NPT 184 
40 Red 8 — 450 
87 eed EEN) A84 

285 47 18 
MEV Re, 86 
47 8 

Rd 0 Wijzergetallen: 

47 55.2, 1, 5,8, 12 NN 

O2 Nee CO 
En 268 
20 sd 
IR ED 
IOR ELO 
174 — 68 
ORDE BELT 
174 LE 
150 4290 
174 300 
TEE Wat 
— 409 — 484 
RT — 409 
—- 893 5484 
— 978 ROS 


950 enz. 


13 


Na het optreden van het bijzonder fortuinlijke wijzergetal *) 490 
kunnen wij de bewerking onmiddellijk staken. 
Wij vinden de naderende breuken: 
on If son Al 196 463: 659” 323373 
Reen GoneslOk 2575 694 -163° 299'« 113843: 
De laatste breuk geeft den grootsten wortel tot op één enkele 
eenheid van de 10° decimaal nauwkeurig: 
Xx, = 2,84051 72035. 
B. Ontwikkeling van den kleinsten positieven wortel. 








2 3 — 18 „4 
— 104 — 16 16 27 
— 102 — 13 — 2 29 
— 220 112 16 208 
— 322 sE 14 237 
— 5634 — 102 16 4985 
— 0956 —d 30 9222 
— 2477 — 102 — 3 — 743 
— 8433 — 105 27 4579 
— 102 — 105 
— 207 — 78 
152 116 
— DO 38 Wijzergetallen: 
616 Oek Oe ende Onhe, 0: 
061 154 
— 644 116 
— 83 270 
— 644 — 166 
— 727 104 
— 644 — 1454 
— 1371 — 1350 
432 1422 
— 939 72 
8964 1422 


en Pa eener: Zie volg. blz. 


1) Het voordeel van het vinden van een zoo hoog wijzergetal springt hier duidelijk 
in het oog. Men losse vooral de gegeven vergelijking ook volgens de methode van 
Newton-Horner op en zal dan spoedig ontdekken, dat de bewerkingen daarbij veel 
meer gereken en zorg vereischen. 

Als een uitkomst in 6 à 7 decimalen voldoende is kan elke gevonden waarde van 
x logarithmisch gecontroleerd worden, omdat steeds 


3 
141  — 263 
eee 





moet zijn. 
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8025 1494 

— 5956 1422 
ie 9060 NOAR 
— 5956 2069 
__ 3887 4985 
20906 0887 
— 9843 1098 
6320 5222 
REDE LANE oi) 
3 11542 5222 
— 9065 11542 
— 8433 5222 
— 17498 16764 
PeRA38 Rt 17408 
OBO EERT Ta R 
BR4A3 Ad Pa 95081 
Ci 5486L u TET 0HBbS 


Uit deze wijzergetallen vindt men als naderende breuken: 














1 gee Bei 191102 490% 1100 VELO 
Ot eet 0, 47e 33e 37 TET ESO 
15298 23497 _ 156280 
5549’ 8523 °“ 56687 


Uit de laatste leiden wij af: x, = 2,756893115. 

C. Den negatieven wortel bepalen wij door in de gegeven 
vergelijking de wortels tegengesteld te nemen. De bewerking ver- 
loopt geheel normaal: 











— 6 | 0 141 263 
eerdre leek pek — 45 
— 187 — 30 —g 218% 
50 tell — 300 104 
— 137 — 60 — 809 322 
— 505197 — 30 218 1170 
— 505334 — 90 — 91 1492 
‚98514 Bel, med L8 170455 
— 603848 — 181 LL 171947 
_17 218 
— 04 345 
Ee meu 
— 241 104 
nie ne 
— 428 — 824 Wijzergetallen: 
— 187 644 


5, 11,2, 15, 2 en 
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— 615 320 
640 644 

95 964 

1928 644 

1953 1608 

— 2055 — 1530 
102 78 

— 055 _-— 32355 
Dot 30977 
— 2055 31332 
— 4212 — 945 
— 19845 31332 
— 24057 30387 
638127 31332 
614070 61719 

— 505334 108736 
108736 170455 

— 505334 « — 396598 
— 396598 _— 226143 
— 505334 515841 
— 901932 289598 
869094 515841 

— 32838 805539 
2416617 515841 


2383779 1321380 
Uit deze wijzergetallen volgen als naderende breuken: 
Bou Ll “28 431 9079 «9510 37609 159946 
Dee 2e Oerd 1622 1699. 6719” 28575" 
De laatste breuk geeft: Xx; = (—) 5,09741032. 








Contrôle. | EPR OD, 
OSL 

X3 = — 9,09741 032 

Oke tekoes e 40, 00000:000 


Tot nu toe hebben wij ondersteld, dat alle wortels der op te 
lossen vergelijking reëel waren. De methode kan echter zonder 
eenige wijziging toegepast worden voor de berekening van een 
reëelen wortel, wanneer de anderen allen of gedeeltelijk complex 
zijn. Ook dan blijft het wijzergetal bij benadering gelijk aan 


oes 


— A Immers, wanneer men alle wortels eener vergelijking ver- 
0 
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mindert met het aantal geheelen van den grootsten positieven 
wortel, zal deze laatste overgaan ín een wortel, die nog positief 
is, maar een getallenwaarde heeft, die tusschen O en 1 ligt; de 
kleinere reëele wortels krijgen een volstrekte waarde, die grooter 
is; wanneer de complexe wortels een geringen modulus hebben, 
zullen zij overgaan ín andere complexe wortels, wier moduli in 
hoofdzaak bepaald worden door het getal, dat er van wordt 
afgetrokken (en hetwelk minstens = 1 is, maar dikwijls veel 
grooter). Bij omkeering der wortels zal de positieve, die tusschen 
O0 en 1 ligt, overgaan in een wortel > 1; de andere reëele wortels 
in dergelijke met volstrekte waarde < 1 (soms zelfs met zeer 
geringe volstrekte waarde); ook de complexe wortels zullen over- 
gaan in wortels van dezelfde soort, maar ín het door ons gedachte 
geval toch weer moduli <1 krijgen en dus evenals de andere, 
reëele, wortels weinig gewicht in de schaal leggen. 


Zoo zal men, na omkeering, weer uit — Ai netuijen kun- 


Ao 
nen opmaken. Wanneer de complexe wortels een grooten modulus 
hebben, zal deze, nadat de wortels met een aantal geheelen zijn 
verminderd, belangrijk blijven, maar na omkeering weer een gering 
bedrag krijgen. Zoo komen wij dan weer op het 1° geval terug. 

Beter dan uitvoerige beschouwingen zal weer een voorbeeld 
de bruikbaarheid der methode kunnen toelichten, ook bij aan- 
wezigheid van complexe wortels. Wij kiezen tot voorbeeld de 
vergelijking: | 

xt — 2x3 — BA — AKAH 5 =O. 

Zij heeft twee positieve wortels, waarvan de grootste tusschen 
3 en 4 en de kleinste tusschen O en 1 blijkt te liggen; boven- 
dien twee complexe wortels, die wij zulien berekenen, als de 
reëele gevonden zijn. 

Berekening van den grootsten positieven wortel. 


1 Ep EA zj 5 
500 8 GORE te alt, 
501 ek 0 B DE 
61872 3 12 vO EN 
62373 MERE BRD 32 — 435 
anne hen — 88 
dj 33 E43 


3 — 15 — 39 Zie volg. blz. 








10 18 18 
Hee ads — 35 
100. — 172 — 73 Wijzergetallen: 
or tlnk #00 VESEEMERENDEN VANG 
— 760 „5 A7 — 108 
Ee A, 106, 
242 — 53 BT 4 
irl a LEE ARE 
1244 2435 4656 
_102 492 __— 5220 
2246 6927 — 564 
_ 1002 — 6768 25290 
3248 159 — 5784 
108 — 6408 __— 520 
5156 — 69249 — 11004 
per le20iE, 
ORL 201207 


Naderende breuken: 
Omnbotnkos dor, 4300 20! 
Klee ee 1 lar 137 E02 
de laatste geeft: x, —= 3,182478. 
Berekening van den kleinsten positieven wortel. 
Omdat het aantal geheelen = 0 ís), kan de Algorithmus hier 
terstond aan de rechterhand beginnen. 























l — 2 — d — 4 5) 
— 4 — 2 l 5 4 
— ò — 4 — 2 1: 49 
— 02 N 6 5) 164 
— 99 0 t 6 213 
— 10460 — 6 Led o 43236 
Ee 10515 — 6 15 Il 43449 
— 1488844 — 6 — 12 5) 
— 1499359 Eje 3 16 
en hm 6 Zie volg. biz. 
1) Dat dit aantal geheelen —= 0 is, blijkt onmiddellijk door het met één vermin- 
deren van de variaties, als de wortels met 1 verminderd worden: 
1 — 2 — 3 — 4 5) 
l — | — 4 — 8 
Eg en IE Ee 
1 0 — + 
ee ed, tl 
1 1 
Pp DE| 
l 
0) 
































AR | 29 
EN; Li35 ED) 
—- 24 BER SI 
es js) 29 49 
ED ed Es 
50 78 49 
Be EN) TE Wijzergetallen : 
BE wel 49 0,152 TREE 
Koast12 AGT schr 
A0 — 924 49 
ia D: UT LTO 
TTO — 444 Ed 
339 — 491 82 
et E10 — 664 — 982 
— 449 — 1155 — 900 
— 1650 825 1065 
— 2092 — 330 165 
29100 6150 1065 
27008 __ 5820 > 1230 
— 31545 11475 1065 
— 4537 17295 2995 
Std ESO 1065 
36082 UR en OGS 3360 
— 31545 _— 108246 11052 
Ls DTO Ta LOD 14412 
31545 — 202881 + == 313686 
00172 807443 209274 
— 1913520 34083 304143 
SC ODI2602 013360 0e RAIO 
13228068 2163084 304143 
11215376 1889724 309012 
4292105 304143 
6181829 613155 
304143 
11917298 


Benaderende breuken: 
l 0 | Dn Bte 040 tT SLEE 
OPP TBA Lod SSN 
Uit de laatste volgt: …x,-=.0,728 1272. 
NaabukREe tT Xe 20,9 LL ZOO EE de 
Met behulp van de gegeven vergelijking vinden wij dan: 


Ker Wen FO 1 305 EEE — 2,155955; 











5 
2,319158 
X3 en X, zullen dus de wortels zijn van: 

Xx. 1,911305x + 2,155955 = 0. 


RD 


Bij oplossing vindt men: 

X3 = — 0,955602 + 1,114800í en x,= — 0,955602 — 1,114800%. 

Er waren dus 2 complexe wortels, maar de oplossing der reëele 
wortels verloopt geheel normaal. Met behulp der gevonden wijzer- 
getallen zou gemakkelijk nagegaan kunnen worden hoe de com- 
plexe wortels in den loop der bewerkingen telkens gewijzigd 
worden en welken invloed zij uitoefenen op de som van alle 
wortels en daardoor op de bepaling van het aantal geheelen, 
waarmee telkens verminderd moet worden }). 

Als verdere toepassing van de methode ter berekening van een 
reëelen wortel bij aanwezigheid van eenige complexe zullen wij 
hier nog behandelen de ontwikkeling van den n°-machtswortel 
uit een gegeven getal A in een kettingbreuk. 

De berekening van x= WÀ kan nl. opgevat worden als de be- 
rekening van den (positieven of negatieven) reëelen wortel van de 
vergelijking x”—ÀÂ=0. Op deze vergelijking passen wij weer afwis- 
selend van links naar rechts en omgekeerd den Algorithmus toe °). 

Voorbeeld. Gevraagd zij 1/82. 


Ì 0 0 — 82 

9 4 16 64 

73 + 16 — 18 

405 ies 92 l 

478 8 48 — 17 

— 297 En — 96 — 1512 

181 12 12 — 1529 
24 — 36 Zie volg. blz 


D Wij behouden ons voor, over dezen invloed later in bijzonderheden te treden. 
2) De gewone Newton-Hornersche methode zou hier geven: 








1 0 0 E39 
4 16 64 
4 16 ER DH 
4 32 15,507 
A 8 48 0 408 
4 3,69 
12 51.69 
0,3 3,78 
12,3 55,47 
0,3 
12,6 
0,3 





12,9 enz. enz. 





36 — 24 
— 48 — 96 
— 12 — 60 
73 61 
61 l 
73 134 Wijzergetallen : 

134 135 4, 2,159 ORO 
78 — 153 
207 — 18 
— 162 — 153 
45 — 171 
— 1539 — 153 
— 1494 — 824 
1434 — 180 
— 60 — 504 
1434 4122 
1374 3618 
3 1434 — 4587 
2808 — 969 
— 2907 — 4587 
— 99 — 5506 
— 16668 — 4587 
— 16757 — 101438 


Naderende breuken: 
Le TAD ALS AEL 200 IO VO 
Oba rleAE DBP 29 ADO LRZ 
Uit de laatste naderende breuk kan men afleiden: 
v82 = 4,34448149 tot en met de laatste decimaal nauwkeurig. 
Bij wvierkantswortels en in het algemeen bij de ontwikkeling van 
wortels van vierkantsvergelijkingen in kettingbreuken kan men 
dikwijls met een zeer beperkt aantal bewerkingen volstaan, omdat, 
zooals bekend is, de te vinden kettingbreuken periodiek (anders 
gezegd repeteerend) moeten zijn. Men staakt de bewerking zoodra 
een vroegere coëfticientenrij is teruggekomen (met dezelfde of 
met juist tegengestelde teekens). 
Voorbeeld. Zij gevraagd y 47 in een kettingbreuk te ontwikkelen. 
Wij zullen de achtereenvolgende coëfficientenrijen naast de 











bewerking opschrijven totdat een vroegere rij terugkeert. Dan 


keeren de vroegere wijzergetallen ook weer terug. 
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1 0 — 47 
l 6 36 Wijzergetal, vol- 
2 6 — 11 Coëfficienten-rijen gend op deze rij, 
— 1 6 Ge vrel, | 
1 12 — 11 l 0 — 47 6 
— 11 0 
| —_ il <-— 11 12 l 1 
— 11 
WD 2 —10 — ll 5 
mest ze 10 2 1 
0 
10 rg ee 12; 
10 
— 11 <- — 11 12 l Ikenz 
Ze 
— 11 
— 12 
12 
0 
E2 
12 


De gevraagde wijzergetallen zijn: 6; 1,5, 1, 12: 
De kettingbreuk is gemengd-repeteerend; de periode der wijzer- 
getallen is 1, 5, 1, 12. 
Voorbeeld. Men vraagt de wortels van de vergelijking: 
Ax3 — 2x +21 =0 
in een kettingbreuk te ontwikkelen. 
1°. Bewerking voor den grootsten wortel. 





4 — 22 21 
Ees BMRO rit Wijzergetal; vol- 
7 — 6 — 3 Coëfficienten-rijen. gend op deze rij. 
— 16 16 — 1 
— 9 10 — 4 4 _— 22 21 4 
A SE 
1 (2) —3 10 4 5) 
mk 
ide 7 — 8 — ò 1 
7 
atd — 4 6 7 2: 
EE, 
6 (5) el 8 enz. 
zee 
— 2 
— 8 
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De coëfficienten-rij (5) is blijkbaar identiek met (2). De periode 
wordtrd nis 12 

De grootste wortel wordt dus voorgesteld door de gemengd- 
repeteerende kettingbreuk: 4; 98, 1, 2. 

20 Bewerking voor den kleinsten wortel. 





4 — 99 21 | 
ie. Nd ke) Wijzergetal, vol- 
— 4 — 18 3 Coëfficienten-rijen. gend op deze rij. 
1 4 4 
ED edet TEA ANNO 21 1 
12 EE 
DEE EEE 4 4 
12 | 
TG REEN Bechet 5 A 
Zn 
WD 13 MENG Seit 1 
EES 
Eer ER 8 7 3; 
d 
Ar. (BR EE) sb 9 enz. 
7 ) 
ET: 
eig 
Dl 
ed, 
AD 


De kleinste wortel kan dus worden voorgesteld door de gemengd- 
repeteerende kettingbreuk: 1, 4; 2, 1,3. 

Beide wortels hebben dezelfde getallen in hun periode; de 
volgorde der wijzergetallen is (hier) echter omgekeerd. 

Als een laatste aanbeveling der methode zij hier vermeld, dat 
men door hare toepassing geleid kan worden tot het opsporen van 
rationeele tweede-machts-deelers van vormen van hoogeren graad. 

Wanneer men nl. bij de berekening van zekeren wortel een rij 
wijzergetallen in dezelfde volgorde ziet terugkeeren, komt men 
tot het vermoeden, dat bij dezen wortel een repeteerende ketting- 
breuk behoort. Men kan dan gemakkelijk de vierkantsvergelijking 
afleiden, waarin deze wortel zou voorkomen. Als het eerste lid 
hiervan inderdaad een factor ís van het eerste lid der gegeven 
vergelijking, heeft men meteen een tweeden wortel daarvan 
gevonden. 
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„Voorbeeld. Zij gegeven de vergelijking: 
2xt — BX3 — IX? + 61 — 59 = 0. 
Voor den grootsten positieven wortel vindt men hier de wijzer- 
getallen: A A A AE 
De coëfficienten-rijen, waaruit deze volgen, mogen hier — ter 
contrôle — worden opgegeven: 


2 ES ed Bis — 55 

lijf: 7 D/ 16 2) 

86 16 fe O0 4) 27 
05 49 369 328 86 
1250 328 RAT — 631 285 

Meal 17 1 631 5139 4672 „01200 


Hoewel geen enkele vroegere coëfficientenrij terugkeert, komen 
wij toch uit het voortdurend terugkeeren van de periode 2, 1 in 
de wijzergetallen tot het vermoeden, dat de grootste wortel voor- 
gesteld zal kunnen worden door de gemengd-periodieke ketting- 
breuk : EE AKI 

De vierkants-vergelijking, waarin deze wortel zou voorkomen, 
kan worden gevonden door eliminatie van x, en X uit onder- 
staande betrekkingen: 


l Ì jl 
(1) EE ae (2) Zere (3) rn 


zat odf 
sarge 
waaruit volgt: (4) ..... Kit — UX, —-2=0. 


X1 


eten A 





(2) geeft dan: x, —2= 





sle” 
Uit (3) volgt: Pan 


Uit (1) leiden wij dan af: x,= 


Door substitutie in (4) blijkt dan: 
Ì 2 
B 
Onderzoekt men nu dezen deeler, dan vindt men terstond de 
ontbinding van het eerste lid onzer vergelijking in een product 
van rationeele tweede-machts-deelers: 
2x* — 8xS — 9x? + 6lx — 55 = (2x? — 10X + 11) (£2 + X — 5). 
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IL. ALGEMEENE BENADERINGSMETHODE VOOR DE REËELE 
EN COMPLEXE WORTELS VAN (NUMERIEKE) 
ALGEBRAÏSCHE EN TRANSCENDENTE 
VERGELIJKINGEN. 


L. Grondbeginsel en Geschiedenis. 


S 1. Om terstond met de deur in huis te vallen: het grond- 
beginsel van de hierna te behandelen methode is de scheiding 
van de wortels eener vergelijking door machtsverheffing. 

Wanneer ok NLD k(k > 1), zal Eis =P; 

Wi, w? 
In het algemeen zal reeds voor beperkte waarden van p de 


grootheid kP? elke willekeurig te stellen grens kunnen overschrijden, 
zoodat w‚? en w‚P, al lagen w‚, en w, dicht bij elkaar, ver uiteen 
zullen liggen en zelfs van verschillende grootte-orde zullen zijn. 

Wij zullen in het vervolg b van hoogere grootte-orde noemen 
dan a, wanneer binnen de eenmaal vastgestelde grenzen eener 
berekening (in 5 decimalen bijv.) a ten opzichte van b kan ver- 


waarloosd worden [des b.v. 4 fe 


Drake OE 
Ze gl 1 1 
Wenscht men nu, dat 7 S or dan moet dus —, S js Of 
5 
5 En 
kes 105, Oe 


In onderstaand lijstje vindt men bij eenige waarden van k de 
vereischte waarden van p (naar boven afgerond) opgegeven. 























cd EDE PE | 14115 J 456 FT, Zoe GO 
p= 214 64} 444135 1729 200016 Zee ZO 
/ mod. ws 
Als w, en w, complexe getallen zijn, zal Hodes k gesteld 
kunnen worden. Dan zal Md vO 





mod. w‚? (mod. w)? 
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De modulus van w‚ep zal dan verdwijnend klein worden ver- 
geleken met dien van w‚? en dus zal w‚? dan verwaarloosd 
kunnen worden t.o.v. van wo?. 


S 2. Op bovenvermelde eenvoudige grondgedachte gelukte 
het mij de algemeene methode op te bouwen, die den titel van 
dit opstel vormt. 

Terwijl ík hare ontwikkeling steeds verder zag voortschrijden 
en haar geldigheidsgebied steeds verder leerde uitstrekken, kwam 
het mij steeds onbegrijpelijker voor, dat eenzelfde gedachtengang 
nog nimmer te voren bij andere wiskundigen zich geopenbaard 
zou hebben. En werkelijk bleek mij bij onderzoek, dat de 
magische kracht der machtsverheffing reeds meermalen met goed 
gevolg ín dienst gesteld was van de scheiding en benadering 
van de (reëele) wortels van vergelijkingen. Zelfs vond ik meer 
dan ík gehoopt had: niet alleen trof ik sporen van dezelfde grond- 
gedachte aan bij Newton !), Daniël Bernouillí ® en Euler ®), maar 
ontdekte ook in Gräffe den gelukkige, die als uitvinder der 
methode mag gelden, — ten minste voor zoover de reëele en 
de moduli van de complexe wortels van algebraïsche verge- 
lijkingen betreft. 

Hoewel Francceur *) als zijn voorganger genoemd moet worden, 
mag Gräffe's naam met recht verbonden blijven aan een methode, 
die geheel beantwoordt aan het ideaal van Lagrange, geformuleerd 
in diens „Mémoire sur la résolution des équations” van 1767: 

„Etant donnée une équation numérique, sans aucune notion 
de la grandeur ní de la nature des racines, en trouver les valeurs 
numériques, exactes sil est possible, ou aussi approchées qu'on 
voudra.” 

De methode heeft geen enkele voorbereiding noodig; de 
theorema's van Descartes, Rolle, Sturm enz. kunnen geheel buiten 
beschouwing blijven. De berekeningen, die zij vereischt, zijn be- 
knopt en van zeer elementairen aard; zij kunnen door eerst- 


1) Arithmetica Universalis T. II, Cap. IV: Sed ad radicem maximam propius accedes 
si quaeras summam quadrato-quadratorum et extrahas ejus radicem quadrato-quadraticam, 
et adhuc magis, si quaeras summam cubo-cuborum et extrahas ejus radicem cubo- 
cubicam. Et ita in infinitum. 

2) Comm. Acad. Petrop III. 

9) Introd. in Analysin Infinitesimorum, Lausannae 1748. 

t) Cours Complet de Mathématiques, $ 58, Paris 1809. 
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beginnende wiskundigen met gemak en haast werktuigelijk verricht 
worden. De aard der wortels (reëel of complex) komt in den 
loop der becijferingen van zelf aan den dag. 

Gräffe zelf ging in zijne verhandeling, gepubliceerd in 1837 en 
in 1839 bekroond door de Berlijnsche Academie, niet verder dan 
de bepaling van de reëele en de moduli van de complexe wortels, 
in de onderstelling, dat deze grootheden onderling verschillend 
zijn. De methode vond spoedig een bewonderaar in den grooten 
astronoom Encke, die haar aanbeval en uitbreidde in een ver- 
handeling, opgenomen in het Berliner Astr. Jahrbuch van 1841 
en later op verzoek der redactie overgedrukt in Band XXII van 
Poggendorf's Annalen (p. 133, s.q.q.). 

Hoewel door dezen herdruk voorkomen was, dat het Berliner 
Ästr. Jahrbuch voor Encke’s verhandeling het vergeetboek werd, 
vond de methode in Duitschland langen tijd niet de waardeering, 
die zij verdient. L. Matthiessen noemt haar alleen in het Litteratur- 
verzeichniss van zijne „Grundzüge der litteralen Gleichungen”'; 
in Weber's „Lehrbuch der Algebra’ wordt zij vluchtig besproken, 
maar daarbij wordt alleen acht geslagen op de reëele wortels van 
algebraïsche vergelijkingen. 

Een Spaansch astronoom, Miguel Merino, kreeg Encke's ver- 
handeling bij toeval ín handen en werd terstond een groot 
vereerder der methode. Zijn vrije vertaling van Encke's werk 
werd met vele wijzigingen en toevoegsels een boekdeel van 
260 bladzijden (1879). 

Na den Spanjaard Merino treedt de Franschman Carvallo als 
geestdriftig partijganger van Gräffe's theorie op. Hij maakt er 
zijn „thèse de doctorat” over, die hij later uitgeeft als: „Méthode 
pratique pour la résolution complète des équations algébriques 
ou transcendantes” (1896). In zijn historisch overzicht spreekt hij 
in de volgende bewoordingen over het beginsel van de scheiding 
der wortels door machtsverheffing: „Ainsi grossies, les racines 
sont séparées et immédiatement mesurables, comme les objets 
fins et rapprochés sont séparés et rendus mesurables par le 
microscope.” lets verder zegt hij: „Gräffe trouve toutes les racines 
par une opération plusieurs fois répétée, comme le micrographe 
augmenterait le grossissement par une série de verres gradués.”’ 

Na ín Spanje en Frankrijk bewonderaars gevonden te hebben 
komt de methode weer tot grootere waardeering in Duitschland: 


VA 


C. Runge heeft haar in zijn „Praxis der Gleichungen” }) een 
belangrijke plaats ingeruimd. 

Na bovengenoemde aanbevelingen zal nu een beknopte uiteen- 
zetting der theorie, naar ik hoop, op genoegzame belangstelling 
kunnen rekenen. 


IL. Met Hoofd-theorema. 


Meen Zij 
Ae Ae dente Anne or An =lOhd Ee (IE) 


een gegeven vergelijking. Wij noemen de wortels, gerangschikt 
naar opklimmende volstrekte waarde of modulus z,, 2, enz... Zn. 
Van (1) zullen wij overgaan op een nieuwe vergelijking, die 
de p° machten van de wortels van (1) tot wortels heeft. 
Dit geschiedt door te stellen: 


oe OT ZP 


Wij zullen later °) bijzonderheden geven over de uitvoering 
dezer substitutie, maar merken terstond op, dat de resulteerende 
vergelijking in y van denzelfden graad zal zijn als de geg. verg. 
mnmers,. blijkens y= 2? (p= geheel getal)-zal. met elken 
wortel van (1) slechts één enkele wortel van de nieuwe verge- 
lijking overeenstemmen, zoodat de laatste wortels in hetzelfde 
aantal aanwezig zullen zijn als de eerste *). 

De resulteerende vergelijking moge voorgesteld worden door: 


Nn eenn Barel Baeoura nr (2) 


en „eindvergelijking’” genoemd worden. 

Als p groot genoeg is, zal elk paar wortels van verschillende 
grootte of van verschillenden modulus van de verg. (1) overgaan 
in een tweetal van verschillende orde van (2), terwijl een groep 
wortels van gelijke grootte of modulus overgaat in een groep 
dergelijke wortels der eindvergelijking. 

Allereerst zal nu worden bewezen het volgende: 


1) Sammlung Schubert XIV 1900. 

2) Zie.S 5. 

3) Bovendien kan men volgens de methode van Sylvester gemakkelijk de resultante 
in determinant-vorm voorstellen. 
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Theorema. Als yrp‚1 van hoogere (grootte-) orde is dan ye, 
valt de eindvergelijking uiteen in: 
vr EE Bjyerbe, CA Bies OE 
en Barn J Baker + Wrs B (DFE AN (4) 
De eerste dezer verg. bevat (afgezien van k wortels = 0) de 
n — k grootste wortels van (2), de tweede de Z kleinste wortels }). 
Bewijs. Als yr‚: van hoogere orde is dan yr, is dit à fortiori 
het geval ten opzichte van ye, Yr—2, «….…. Yr 
Nu is: 
Bi Yn ARA AT Pe Ve dT Ve Het TE +}; 
Bi Vn Vitis okeee Mae Meseta ]; 
— Bi == Vn Ira Varian dn Wera Vree Ve etna } 


me ne a ee a en en ee ne 


— IB = Yn Inl «Yarr + | ]. 

Tusschen [ ] denke men telkens vereenigd de termen, die 
blijkens de praemisse van het theorema van lagere orde zijn. 
Doet men dit, dan blijken terstond yn, yn—i, . ….. Yr+1 de wortels 
te zullen zijn van 


pek AAD VE Tr ne RO 
Vit Bidden liae rn chaDERnE 
Dat verder | 

Bie vt Bassant seen ek 
juist de wortels yr, yr, .... #, bevat, kan bijvoorbeeld blijken 
door (2) op omgekeerde wortels te nemen. 


Noemen wij 7 n, dan ontstaat 


Brinke B sn Ae desta Metel VR 








Door, de omkeering zijn. nu” #,, Vas. VS Vern 
gegaan in Ans Dn—1ls olie tels Dn—=k+1s EE BON 
jl ] 
OBN IS ng a el Mi 
Vk Vk+1 
k AE 
Als vroeger ERE Ö,. zal nu == den dUss nas 
Yk+1 Dn—k+1 


hoogere orde zijn dan np. 





1) Een wortel zal hier steeds grooter heeten dan een andere, als zijn volstrekte 
waarde of zijn modulus dien van den anderen overtreft. 
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Volgens het reeds bewezen eerste deel van het theorema zal 

dan de vergelijking: 
Ben + Brunt +... + Brant” = 0 
BNS WOLLGIS 1,5 Agis oee An=kar DEVALEN: 
Door nu weer om te keeren vindt men dat 
Binn ke0 A5 Bransen fr HER zj B =d, 
iist de wortels y,, ys, .... Ve bezit. 

Hiermede is het gestelde in zijn geheel bewezen. 

Opmerking. In het bewijs wordt alleen ondersteld, dat yr. 
vansshoogere orde is dan ye; onder vj, .……..9r EN Yrais Veeas 
……..Yn mag een willekeurig aantal groepen wortels van gelijken 
of vergelijkbaren modulus voorkomen. 


S 4. Beschouwen we nu eens een groep van / wortels, volgende 
OP Yr, TUS Yarr, Yra2, ««..Yrer 

Onderstellen wij yr.1 van hoogere orde dan yr, maar y#+: van 
lagere orde dan yrs141- 

Uit het theorema volgt dan: 

1°. en eed ar Bm Vet ERN ols Pld. 
heeft tot wortels: 


VREDEN ALSI, stare VR VRerin sees Vis 
B een Ven et 0 
heeft de wortels: | 
| Nl VRSle ones ce Vis 
Een E WOIEls Weers Vaer tue. Verl 


zullen juist bevat zijn in: 
Ne henk ek arn tn Bas V ij, 

Derhalve: 

Voor een groep van l wortels, waaraan een wortel van lagere 
orde voorafgaat en waarop een wortel van hoogere orde volgt, 
kan dus een groep van l—+ l opeenvolgende termen uit de eind- 
vergelijking gelicht worden. 

Door deze eigenschap kan men de wortels van gelijken modulus 
afzonderen. Zij is echter ook van groot nut in het geval, dat er 
in een vergelijking een beperkt aantal wortels van weinig ver- 
schillenden modulus voorkomt. Deze wortels zullen nl. slechts bij 
verheffing tot een macht met belangrijken exponent p tot ver- 
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schillende orden gebracht kunnen worden, terwijl dit voor de 
overige wortels wellicht reeds voor een zeer bescheiden waarde 
van p mogelijk zal zijn. Men kan zich met behulp van deze 
eigenschap dus dikwijls zeer beperken bij de keuze van den 
exponent p. 

Het bovenstaande wordt toegelicht in voorbeeld IL. Eerst dienen 
echter enkele bijzondere gevallen van het theorema te worden 
beschouwd en de middelen besproken te worden, waardoor de 
eindvergelijking op de eenvoudigste wijze wordt gevonden. 

Hebben alle wortels verschillenden modulus, dan valt de eind- 
vergelijking uiteen in: 

ya Br= 0, y= Be. 0, Bi bet 0 Nn 
B-2y nn les —= 0) en By ar B, = ON 


De kleinste wortel is dan py, = — ee de grootste y„ = — B; 
n—l 
in het algemeen zal 
Bn + 
Vk EL EEE ee e je k 5 7 5 = (5) 
Uit ye zal dan volgen: 
P ej 


hr Marijne 
Hieruit vindt men dan onmiddellijk den modulus van ze; het 
argument vereischt nader onderzoek. Immers, als het argument 
van Yr, afgeleid uit (5), = © is, zal dat van ze = es 2e, als 
m een der waarden 0, 1, 2, ....p—l BRK 
Het argument van ze blijft dus voorloopig tot op eenige malen 


ns 
— onzeker. 
P 


De middelen ter nadere vaststelling van dit argument zullen 
later worden besproken. 
Zijn alle wortels twee aan twee deconiugeeratonn (of gelijk), 
dan zullen y„ en y„—: volgen uit: 
By + B, = 0; 
verder vaas en. yri-s ite 
By? + Bsy + B, = 0; 
en zoo vervolgens; 
de beide kleinste wortels y, en yo zullen bevat zijn in de verg.: 


Biz y? jr By in end NO: 


Call No. 
DUPLICATE 


he: \ 
gk ed ed ves TA 
/ ze AJ ep ese veed WE bad eneen ser dopaehi art nT heteen hab dn Amis db idee re audlwes vosvecee Pe 


Methoden on het 


aische en transcende: 


Groninsen P.-N vordhoff 
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UL. Bijzonderheden. 


S 5. Het uitvoeren der transformatie. 
Door te substitueeren z = v/z! leide men uit (1) de vergelijking 
Star Har d.d Han=0. . (6) 
af, welker wortels de quadraten zijn van die der eerste. Deze 
vergelijking kan de „eerste getransformeerde” genoemd worden. 

Daarna levert de substitutie z/= Vz’ een tweede getrans- 
formeerde, die de 4° machten der wortels van (1) tot oplossingen 
heeft. Op deze wijze kan men steeds verder gaan. Na 5 transfor- 
maties is er reeds een vergelijking ontstaan met wortels, gelijk 
aan de 32° machten van die der oorspronkelijke vergelijking. 

Reeds Gräffe heeft gewezen op de eenvoudige vormingswijze 
van (6) uit (1). Men kan de rationeele termen in het eerste lid 
vereenigen, alle termen, die een factor 2’ bevatten, in het tweede 
lid vereenigen, daarna beide leden quadrateeren en eindelijk de 
termen met gelijknamige machten van z’ vereenigen. 

Ik zou hier echter, — vooral met het oog op latere toepassing 
en uitbreiding, — de volgende afleiding willen voorslaan: }) 

Laat (1) voorgesteld worden door F(z) = 0. 

De vergelijking F(2) X F(—z)=0 heeft behalve de wortels 
van F(z)=0 nog evenveel andere, juist gelijk aan de tegen- 
gestelden van die van (Il). Zij moet dus voorgesteld kunnen 
worden door o(z?) = 0. 

Vervangt men nu z? door 2’, dan ontstaat o(z’) = 0 en deze 
vergelijking moet blijkbaar (6) opleveren. 

Bij vermenigvuldiging van: 

PousAret ln ie, ve irtAn reeet nar Är== 0 
met 
ee Agetndebe nt (IA rele vet (rl )An Zr Än=0 
vindt men voor den coëfficient van z** vooreerst: 
(— IFA? nl. uit het product der beide termen met z*; daarna: 


WTA iAn-rei Uit (— 1) tlAr pikt IX Arpeizk tk; verder: 
En Arrr inr Uit (— ITA eri Th X Aneri2ftL 
De beide laatste termen kunnen vereenigd worden tot: 
Nl 1 Ar priAnekr1- 


1) Het is mij gebleken, dat deze afleiding ook bij Runge te vinden is. 
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Men ziet nu terstond in, dat de coefficient van 2 in F(2) XF(—2)=0, 
dus de coëfficient a, van z’* in (6) zal voorgesteld worden door: 
(Ari 2X (IA inkt 2 X(— 1) Areno 
of door: 

(EHA re VA reine ARR en 

Opmerking. Als men in (6) z/ vervangt door — 2/, zal deze 
coëfficient juist worden : 

(UMA ZA An ett ZA Ar 
of | 
Ark SQA Antas DAA on 

Omdat door de machtsverheffing het teeken, en in het alge- 
meen het argument van de wortels toch in zekeren zin onbepaald 
wordt, verdient het aanbeveling, in het vervolg onder de „eerste 
getransformeerde’ te verstaan de vergelijking, die de tegenge- 
stelden van de quadraten der wortels van (1) tot oplossingen 
heeft. Voor de berekeningen van hare coëfficienten geldt dan 
de volgende: | | 

Regel. Een willekeurige coëfficient van de eerste getransfor- 
meerde is gelijk aan het quadraat van den overeenkomstigen 
coëfficient der oorspronkelijke vergelijking, min het dubbele 
product der beide coëfficienten, die den eersten omsluiten, plus 
het dubbele product der coëfficienten, die de beide vorige insluiten 
enz. totdat men een der beide uiterste termen der egen ver- 
gelijking bereikt heeft. 

De coëfficienten der verdere getransformeerden vindt men dan 
door steeds denzelfden regel op de coëfficienten der vorige ge- 
transformeerde toe te passen. 

Men kan dezen regel gemakkelijk in een beknopten alge- 
braïschen vorm brengen. 

Bij OE At 

k=0 \ 
behooren als coëfficienten der eerste getransformeerde: 
m=n—k 
Ak = heben 1 Arm Aram 
(A's met negatieve indices = 0 te rekenen); 
de coëfficienten der volgende getransformeerde zijn. dan: 


m=n—-k 
m=k—n 


„38 


_S$6. Het eindigen der transformaties. 

Vergelijking van de laatst gevonden getransformeerde met de 
voorlaatste wijst onmiddellijk uit of de eindvergelijking reeds 
gevonden is. Zij nl: | 


Bt Birt Bip t eer Bary Ba 0 
een eindvergelijking, waarbij y‚+1: van hoogere orde is dan y,, 
dan zal deze uiteenvallen in: 


Veel ese Bao Ven 01 ED 
BM “E Bier + ME Ie + B == 0), 
Hierin is 
(— Bao sr Vn Vri seien Vi+1 an [ |. 
Wordt deze eindvergelijking nog eens getransformeerd door 
y= Vv. tot 
VECHT + TEN ene ee Cone ST (A 
dan zal in (7) ‚1 van hoogere orde zijn dan y’‚ en (7) dus 
uiteengenomen kunnen worden in: 


ks, + oro ee Ve == (en: 


ade Ae Ci =D. 
Nu is 
BN en or Mam Merdelete ss M perd tr | | of 
OV at lan Vetere}. 
Men ziet dus, dat 
Gael ol 


De eindvergelijking is dus gevonden; wanneer bij nog een- 
maal voortgezette transformatie een zoo groot mogelijk aantal 
coëfficienten in hunne guadraten overgaan. 

Zijn alle wortels ongelijk, dan gaan alle coëfficienten in hunne 
quadraten over; zijn ze twee aan twee gelijk of geconjugeerd, 
dan ondergaan de coëfficienten deze verandering juist om den 
anderen enz. 

Met het oog op de overwegingen in $ 4 kan men dikwijls in 
plaats van met een zoo groot mogelijk aantal volstaan met een 
genoegzaam groot aantal coëfficienten, dat bij nog eenmaal voort- 
gezette transformatie gequadrateerd wordt. 

Het aantal transformaties. Een groot voordeel der methode is, 
dat dit aantal in het algemeen zeer beperkt kan zijn. Na vijf 

3 
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transformaties zijn alle wortels reeds tot de 2%-macht, d.í. tot de 
32°*-macht verheven; na 8& transformaties reeds tot de 256°- 
macht enz. 

Ter illustratie zullen eenige voorbeelden volgen. De rekening 
geschiedt in 5 decimalen; de (moduli van de) wortels zullen dan 
met ongeveer even groote nauwkeurigheid bekend worden. 

Immers, al gaat bij de verheffing van de wortels tot een hooge 
(bijv. 64°) macht door de weglating van decimalen de absolute 
nauwkeurigheid der coëfficienten verloren, door de latere wortel- 
trekking wordt de zoo ontstane fout (nl. de relatieve fout) weer 
in dezelfde mate verkleind. 

Wanneer men bijvoorbeeld Ri de berekening van 


alit Th Ee 105 Fars te) 














Trekt men hieruit weer den 64°-machtswortel, dan vindt 
men terug: 





64 32 32164 1 32 63 
zij/| K nt a zal! jn ml allt ige zo tes 


De fout, ontstaan door de machtsverheffing, wordt dus te niet 
gedaan door de worteltrekking. 

Het ís verwonderlijk, dat nergens nadele gewezen wordt 
op den gunstigen invloed van deze worteltrekking op de relatieve 
fout van de eind-uitkomst. Encke zegt zelfs met eenigen nadruk, 
dat men wegens de onnauwkeurigheden, aan de transformatie 
verbonden, nooit moet verwachten een gewenschte nauwkeurigheid 
te kunnen bereiken door de berekeningen in het daaraan beant- 
woordend aantal decimalen te verrichten. Eén wenk schijnt mij 
echter behartigenswaard nl. deze, dat men bij de eerste of de 
twee eerste transformaties in sommige gevallen met 7 decimalen 
moet rekenen, als men de einduitkomst in 5 decimalen nauw- 
keurig wil hebben. Door de aftrekkingen Kunnen nl. wel eens 
enkele decimalen verloren gaan, waardoor de verdere berekeningen 
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minder nauwkeurig zouden worden. Maar ín het algemeen is 
deze voorzorg overbodig. Het uit het oog verliezen van den 
gunstigen invloed der finale worteltrekking is de oorzaak, dat 
Encke en ook Carvallo meermalen een uitkomst vinden met 
belangrijk grootere nauwkeurigheid „dan te verwachten was” }). 

Dit laatste mag natuurlijk alleen bij uitzondering voorkomen, — 
als men nl. bij het stellen van zijne verwachting werkelijk alle 
omstandigheden, die fouten kunnen veroorzaken, behoorlijk in 
rekening gebracht heeft. 


Voorbeeld 1: 
ze — J.125 +-0.5722 + 0.1942 — 0.042 = 0. 


De coëfficienten der achtereenvolgende getransformeerden zijn: 


l 9.27 1.0557 0.085516 0.001764 
1 nl. 093357 nl. —0.25278 nl. — 2.44507 nl. — 5.05700 
1 „ 1.86050 „ —0.60097 „ — 5.026838 „ — 11.01400 
l „ 3.72096., —1.21147 „ — 1007819 „ — 22.02800 
1 „ 7.44192 „ — 242305 „ — 2015710 „ — 44.05600 
l „ 14.88384 „ — 484610 „ — 4031420 „ — 88.11200 


Aan de logarithmen ziet men terstond, dat de coëfficienten der 
laatste vergelijking allen gelijk zijn aan de quadraten van de 
overeenkomstige coëfficienten der voorlaatste. Elk van beide kan 
dus als eindvergelijking gelden. De laatste valt uiteen in 4 ge- 
deeltelijke vergelijkingen, waaruit men vindt: 


— 88.11200 + 40.31420 — 025316 — 1: 
64 
hieruit volgt mod. z, = 0.17913. 
Verder blijkt 


mod. z, = 0.27913; mod. 2; = 0.49172 en mod. z, = 1.7083. 


De moduli zijn alle verschillend; daar de coëfficienten der 
gegeven vergelijking alle reëel zijn, kunnen er geen complexe 
wortels zijn. 

Door te letten op de som der wortels (+ 2,1) ziet men onmid- 
dellijk, dat de fweede wortel negatief moet zijn; de andere zijn 
positief. Men heeft dus: 


B RAILS, Zo — 0.219135, 2, = 049172, Enz, = 1.7083. 


log. mod. z, = 


1) Méthode pratique etc. p. 20: 
_ „On le voit, le calcul se trouve plus précis que nous n'avions demandé’ enz. 
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Contrôle. De vergelijking werd gevormd uit: 

(z? + 0.1z — 0.05) (2? — 2.22 + 0.84) = 0. 

Door oplossing van beide, hierin bevatte, vierkantsvergelijkingen 
kan men zich terstond overtuigen, dat de boven gevonden wortels 
allen tot in de laatste decimaal nauwkeurig zijn. 

Voorbeeld II: 

23 32 20. 

De transformaties worden hier voorgesteld door het volgende 

schema: 


zb! 2 z 
DE edt 18 — 8 +4 
De werd + 185 — 40 16 
25e lever tlads0803b Meen 903040 Ted 
Pin 23) „907072. > trip 00. 04048 200 ER 


Er is nu een genoegzaam groot aantal coëfficienten der laatste 
vergelijking gelijk aan het quadraat van een overeenkomstigen 
coëfficient der voorlaatste; de eindvergelijking is gevonden. Zij 
valt uiteen in: | 

y + num. log. 9.07072 = 0 en 
num. log. 9.07072y? + num. log. 6.04048y + 256° = 0 .... (8). 


Deze laatste vergelijking heeft blijkbaar twee complexe wortels; 
de oorspronkelijke dus ook. Voor hun gemeenschappelijken 
modulus vind men: 

mod. 21,2 = Ris = 0.73630. 

Voor enkele lezers is misschien hier eenige toelichting ge- 
wenscht. Als de vierkantsvergelijking az? + bz + c = 0 twee com- 
plexe wortels z, =p + gi en Z, =p — gi heeft, is hun gemeen= 


schappelijke modulus Rio = V(p° +g°®) = VZ = Von | 





EEA 256° 
Daarom is in (8) mod. yo = fee log. 907072: 
32 
RA? 
Hieruit volgt dan mod. zi2 = Ae 


_ 2log. 256 — 9.07072 
ig Dz 
en derhalve Rio = mod. zis = 0,73630. 


dus log. mod. zz — 0.86705 — 1 


1 De machten van 2, die vooraan staan, geven aan tot welke macht de wortels 
der oorspr. verg. verheven zijn. 
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Verder zal y; = num. log. 9.07072 (afgezien van het teeken) en 


16 
AUS. Z; = V num. log. 9.07072; 
IOS AL —= 0,56692; 
16 
teres 


Met één blik op de gegeven vergelijking is het duidelijk, dat 
het teeken van z; positief zal zijn. - 

Het argument der beide complexe wortels zou Kier kunnen 
volgen uit: Ais COS: Ger d00R leer ec Of: 50 safe ve rekt (9) 

Ri COS. Er ar 0.34455, 
hetgeen oplevert: 
| Pel ADW Den los == 242 DEM 

Toelichting. Als 
| Zie =P £ qi = Rip(cos.p + isin. 0), 

is Zin Ze RI COS P. 

Verder moet z, +2, +23 = 3. 

Hieruit volgt terstond (9). 

Opmerking. De eerste twee transformaties zijn rechtstreeks met 
de coëfficienten berekend. De verdere berekeningen zijn verricht 
met gewone logarithmen in 5 decimalen. Men kan ook met voor- 
deel Gaussische logarithmen gebruiken (o.a. Runge); bij be- 
rekeningen, die niet verder dan 3 decimalen behoeven te gaan, 
kan de rekenliniaal goede diensten bewijzen (Carvallo). 

ed UT: 

_— 13:125- 433. 4z? — 12642 + 960 = 0. 

Deze BE. is ontstaan uit: 

(2? — 3.12 + 2.4) (2? — 102 + 400) = 0. 

ERRREeIte daarom twee reëele wortels “zie 1,9 én: 201,65 
bovendien twee complexe, waarvan R34 = 20. 


Na 3 transformaties, — zoodat z = Vy, — zullen de waarden 
van 1.9® en 1.6% nog vrij dicht bij elkaar liggen, maar beide 
reeds van lagere orde zijn dan de andere wortels y;, en ys, 
Reeds na 3 transformaties moet men daarom de vergelijking ín 
twee vierkantsvergelijkingen uiteen kunnen nemen. Dit blijkt aldus: 


pk ed eld z 
Pi l — 2.84210 5.19490 + 5.88400 5.86454 
28 l P5. 20008 10.40828 11,47337 11.92908 
Jh l — 10.34793 20.81649 22.65264 23.85816 


38 


N.B. In plaats van de coëfficienten hebben wij hier overal 
hun logarithmen genomen, voorzien van het teeken, dat bij den 
coëfficient behoort; — 10.34793 beteekent dus — num. log. 10.34793 
enz. De eindvergelijking valt uiteen in: 

y° — 10.34793y + 20.81649 = 0 en 
20.81649y? + 22.65264y + 23.85816 —= 0. 

De eerste hiervan heeft blijkbaar 2 complexe wortels. De 
discriminant is n.l. num. log. 20.69586 — 4 X num. log. 20.81649 = 
num. log. 20.69586 — num. log. 21.41855, dus negatief. 

Men vindt: 


Ran A log. 20.81649, waaruit 
log. R34 = 1.30103 en dus R34 = 20.0000. 
De tweede vergelijking of 
y? + 1.83615y + 3.04167 = O (nog steeds num. log.!) . .. . (10) 
kan nu afzonderlijk verder behandeld worden. 
Hier kan zij echter wegens haar lagen graad rechtstreeks 
worden opgelost. Men vindt: 
Vi =— 25.6350; Yo = — 42.9380. 
Van deze teekens kan men nu afzien. Men vindt dan verder: 


8 
z, = V25.6350 = 1.50003 en 


dn V42.9380 == 1559994 
Door substitutie van z, = 1.5 en 2, = 1.6 in de oorspronkelijke 
vergelijking overtuigt men zich hank dat het teeken van 
z, en 2, positief is. 
Wanneer men op (10) de transformatie verder voortzet, heeft men: 
Ze jl 3.39808 6.08334 
En | 6.58330 12.16668 
dh 1 13.06967 24.33336 
Sd 1 26.12555 __48.66672 
ZE Ì 92.25089 97.33344 
Men behoeft niet verder te gaan. Na de eerst-volgende transfor- 
matie zullen alle coëfficienten overgegaan zijn in hunne quadraten. 
97.33344 — 52.25089 








Men vindt nu log. z, = — 556 =017, 0408 
Zie 1000008, 
52.25089 
Verder loen == OEE 0.20410; 


zo = 1.59904. 
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Opmerking. Bij het quadrateeren van complexe wortels zullen 
hunne argumenten steeds verdubbeld worden en dus afwisselend 
in allerlei quadranten terecht kunnen komen. Op de aanwezigheid 
_ van (geconjugeerd) complexe wortels zal men dus vooral opmerk- 
zaam worden gemaakt, wanneer zich bij sommige coëfficienten 
in den loop der bewerkingen teekenveranderingen voordoen. In 
de beide laatste voorbeelden valt dit gemakkelijk te controleeren. 


IV. Het argument der complexe wortels. 


S 7. Encke leidt uit de eindvergelijking alleen den modulus R 
van zekeren wortel der oorspronkelijke vergelijking af. Het 
argument @ kan dan bepaald worden door substitutie van 
R(cos. 9 + isin.) in de oorspronkelijke vergelijking. Uit deze 
volgen dan, ook bij complexe coëfficienten, twee vergelijkingen 
met @ als onbekende. Omdat cos. ko en a geheel in machten 
van COS. kunnen worden uitgedrukt, zijn deze twee vergelijkingen 
als algebraïsche vergelijkingen in cos. op te vatten, wier coêifi- 
cienten bovendien reëel zijn. 

De grootste gemeene deeler der eerste leden, gelijk nul gesteld, 
levert de waarden van cos. ©, die bij den beschouwden modulus 
R behooren. Is deze modulus enkelvoudig, dan moet de G.G.D. 
van den len graad zijn. Bij aanwezigheid van eenige wortels met 
gelijken R, maar verschillend argument, is deze G.G.D. van 
hoogeren graad. Deze methode zal hier eenigszins nader worden 
toegelicht bij een algemeene vergelijking met reëele coëfficienten, 
vooreerst omdat men in de praktijk alleen met zoodanige te 
doen heeft en ten tweede omdat daarbij zeer belangrijke vereen- 
voudigingen optreden. 

Stelt men in de vergelijking: 


AA enke tnt dese tr Aon Aon 0 
z= R(cos. @ +risin.g), 
dan moeten de bij R behoorende waarden van o tegelijkertijd 
voldoen aan: 
AQRcos.2notA, R?r-icos.(2n—1)pt….tAon—iReOs.otA2n=0.…… (11) 
en 
AoR*sin. 270 + A,R?r-tsin. (2n—1)p +. HAon—tRsin.p =0..… (12) 
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Nu zal cos. ne X (11) + sin. ne X (12) opleveren: 
AoR? cos. no + A,R2-1 cos. (n — Ig +... + 
+ AoniR cos. (n — 1e + Aoncos.ng =0 ....(13) 
Verder geeft — sin. no X (11) + cos. no X (12): 
AoR?-sin. Ao A, R®Esin. (7, —41)G Ste tinten 
— Aon—aiRsin. (n — 1)e — Aonsin.no = 0 .... (14) 
In (18) en (14) kan men de symmetrisch ten opzichte van het 


midden geplaatste termen vereenigen. 
Noemt men: 


AR ie Aon = Bo, AR” a: Aon ges Kh 
A,R-1 + AoniR=B,  AjRI — AopiR = Co 
AR? + Aon—oR* =B3,  AgR? — AonoR? = C5 


mmm ne ee 


AR” = B, An—iR"*1 — Ani?! = Can 
dan gaan (13) en (14) over in: 
B‚cos.ng + B,cos.(n — 1)p +... + Br-1C0S.9 + B =O en .… . (15) 
Cosin. no HG, sin. (n — 1)p 4. + Cor Sin. p=0 (16) 


Om nu tot den enkelen hoek over te gaan heeft men de 
formules noodig: 


n n(n—3 
COS.ng = a OEE — To cost gee 2n-5COS" Gn 








nn — A(n — 5) n—7 n—6, n(n— o)(n— 6 — 7) n—9 ng 
OENE 
en 
En 2rt cos. "lp — ml COS: Aan 
ze En Ne Dn ab: COS, 250 — Dr) erosie 
RER eee 3 
nt Lhee e (nt 00E DRO oo ern 
EVRM RL CAN TET i 29 COS. A ses OS 


De ontwikkelingen (17) en (18) moeten zoover worden voort- 
gezet tot men op negatieve machten van cos.9 zou stuiten. 

Uit (15), (16), (17) en (18), — in de beide laatste betrekkingen 
moet natuurlijk # achtereenvolgens door n — 1, n — 2 enz. ver- 
vangen worden, — kan men nu een resultaat afleiden van de 
gedaante: 
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Bocos.ro HB, cos." IgH B, COS"? HBr COS. 0 Hb =O. (19) 
MOO OS: Ere, nt Hei COS P rr Yar 07 7, (20) 

Bij een vergelijking van den 2nden graad behooren dus twee 
argument-vergelijkingen respectievelijk van den nen en den „ — len 
graad. Voor eenige speciale waarden van 2n tot en met 2n = 8 
worden hieronder de noodige opgaven gedaan. 

Bij vergelijkingen van onevenen, bijv. van den 2n — len Be 
blijven alle bovengenoemde herleidingen — mutatis mutandis — 
gelden. Men krijgt als resultaat weer twee vergelijkingen ín cos. 0, 
die evenals (19) en (20) van den ren en den n — len graad zijn. 
In plaats van dít gewijzigde resultaat afzonderlijk af te leiden 
verdient het de voorkeur, de vergelijking van den 27 — len graad 
door vermenigvuldiging met z tot eene van den 2ren graad 
herleid te denken. 


S 8. Bijzondere gevallen. 
Bij een vergelijking van den 4en graad heeft men de formules: 


AOR AD, AoR* — A, = Co, - 
ere NE zi AARS ERNA Ci 
AakeeiB 
OCOS BEOS Datar B 0 tsern (21) 
DONdOSRe CHE OP ee PA (22) 


De laatste betrekking is op zich zelf reeds voldoende. 


AIEE, 
EOC ee ARS AR 


Hieruit volgt voor een verg. van den Zen graad: 





COS: ee np kt ARE . an een Ue (24) 


Met behulp van (24) vindt men bij het boven behandelde 
voorbeeld II: 
DE D  AMOL ENDE sL). 
Deze uitkomst verschilt 1,1 met de vroeger gevondene; deze 
kan echter, omdat er toevallig voor 2Ri» cos. « slechts 4 decimalen 
overbleven, wel tot op meer dán 1 sec. fout geweest zijn. 
Bij een vergelijking van den 6en graad gelden de formules: 


ARS As —= Bo, AoR* — As = Co, 
AR? ei AsR oer B,, AR? Ae AR a Gj, 
A De, A9R* — A/R* = Ca. 


A;R* = B 


A nnn nn nd 
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4B, cos.*9 + 2B, cos.” @ — (3B, — B) cos. 9 — B, + B; =0; . . «. (25) 


AG, cos.*9 + 2G, cos. ® — GC + C‚ = 0. ee 
Bij die van den 8en graad moet men aanwenden: 

AoR® + As = Bo, AR? — Às mr Er 

ACR Ere AR? — A,R = G,, 

ARP TTAER Sr ARS — AGR? = C 

ASR® TWAS A;R° — A;Rè = 

ARB 

8B, cos.*v + 4B, cos. o — (8B, — 2Bj) cos? # — 
— (3B, — Ei cos. @.-°Bo— Bj Bj = ON (27) 


8C, cos. v + 4C, cos.” 9 — (4C, — 2C9) cos. 9 — CG, + Cz = 0. 
Voorbeeld Iv. En | 
x8 — 9,6x° + 49,92x* — 150.128x% + 320.0528x° — 
412.10176x + 314.987904 = 0 5). 


Op deze verg, welke drie stellen geconjugeerd-complexe wortels 
bezit, — hetgeen bij onze behandeling van zelf aan den dag 
komt, — moge nu de volledige methode toegepast worden. 

De (logarithmen van de) coëff. der ste getransformeerde zijn 
met behulp van 7-decimalige logarithmen bepaald en daarna 
weer bekort tot 5 decimalen. De En verdere rekening ge- 
schiedde ín 5 decimalen. 


— 9.6 49.92 — 150.128 320.0528 — 412.10176 314.987904 
— 088536 239734 — 332905 ‘400630 — 45023000 
— 264378 469768 — 5.35245 7.22575 — 9.00099 999318 


+ 4.97399 936559 — 11.88141 1491010 —+17.82832 19.98636 
+ 9.92635 18.74272 — 24.48679 30.32957 35.47129 39.97272 
+ 18.83467 37.48580 — 49.15259 60.81168 70.67669 79.94544 
— 37.16163 74.97160 — 98.28975 121.63715 141.04648 159.89088 
149.94320 — 196.63616 243.27440 281.75134 319.78176 


Na 7 transformaties is de eindvergelijking gevonden. 

De le, 3e, 5e en 7e coëfficienten zullen bij een volgende trans- 
formatie allen ín hunne quadraten overgaan. De 3e coëif. doet 
dit reeds na de 5e transformatie; daarna kan dan ook reeds een 
drie-termig deel uit de vergelijking gelicht worden. Voor de 
overige coëff. moet men nog enkele transformaties doen. De 
eindvergelijking valt dan uiteen in: 


1) Dit voorbeeld wordt volgens een geheel andere methode behandeld door 
Dr. A. Kempe in Jaarg. XI van het Wiskundig Tijdschrift. 
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(119) py? — nl. 3716163 X y + n.l. 7497160 = 0; 
(29) _n.l. 149.94320 X y? — nl. 196.63616 Xy + n.l. 34327440 =O; 
(39) _n.l. 24327440 Xp? + nl. 281:75134 Xy + n.l. 31978176 — 0. 


Het verloop van den Zen, 4en en Gen coëfficient gedurende de 
transformaties wijst duidelijk op den complexen aard van alle 
wortels. Door de eerste transformatie krijgt men trouwens reeds 
de guadraten van de wortels der oorspr. verg. Als deze dus 
reëele wortels had, zal de le getransformeerde daarvoor alleen 
positieve wortels krijgen. Deze worden dan alleen tegengesteld 
genomen, zoodat de le getransformeerde verg. alleen negatieve 
reëele wortels kan hebben en dus alleen factoren y + y,, y + ya 
enz., als y, en y, positieve getallen voorstellen. Door de reëele 
wortels der oorspronkelijke vergelijking kunnen er dus alleen 
positieve teekens bij de coëfficienten optreden; negatieve teekens 
wijzen op complexe wortels. Ook het feit, dat sommige coëf- 
ficienten geen neiging vertoonen tot quadrateeren en zelfs plot- 
seling van teeken veranderen, wijst op zulke wortels. 

Maar zekerheid geeft hier de eindvergelijking. 

In de gedeeltelijke eindvergelijkingen blijkt de discriminant 
terstond negatief te zijn. | 


Bij (1°) is log. p? = 2 X 37.16163 = 74.32326, 
log. 4g = log. 4 + log. q = 0.60206 + 74.97160 = 75.57366, 

dus p“ — 4q < 0. 
BINZ NIS "log. b° = 2:,196,63616 —= 393.27232 

log. 4ac = log. 4 + log. a + log. c —= 393.81966, 

dus b? — 4ac < 0. 

Bij (3°) is log. b? = 563.50268, 
log. 4ac = 563.65822, 
dus weer b? — 4ac < 0. 

De modulus van beide wortels van (1°) is gelijk aan den wortel 
74.97160 
mins 
Omdat (1°) ontstaan is uit de oorspr. verg. door de wortels tot 
de 64e macht te verheffen (2%), heeft men voor den modulus van 
twee wortels der oorspr. verg. 


log Rs = ne — 0.58571 56, 


uit hun product; de log. van dien modulus is dus 
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Verder heeft men: 
243.27440 — 149.94320 








JOS 556 — = 0.36457 50; 
log. R‚ = 319.78176 — 243.27440 —= 0.29885 6. I) 
‚256 
Voor de berekening van het argument bij R heeft men hier: 
AcR4° = 3268.08 A,R;S = — 81442 
Arte 314.00 AR; = — 1587.52 





B, =3583.07; C‚= 2953.09. B, = — 97317; C,=-—6556.7. 
AR: 10993,3 | 
A‚R;’ = 47496 | 

B, = 157429; C, = 6243.7; «B, = A;R‚3 = — 85824 

Na deeling door den len coëff. (4B, en 4C)) krijgt men de 
argument-vergelijkingen: 

cos. @ — 1.35800 cos.* @ + 0348417 cos. v + 0.080188 — 0 

en 
cos. 9 — 1.11018 cos. # + 0.278573.= 0. 
De rest der deeling van de laatste verg. op de voorlaatste, 
gelijk nul gesteld, geeft: 
— 0.205281 cos. ® + 0.149225 = 0, waaruit: 


0.149225 
0205281’ log. cos. p —= 9.86149; 


© = 43°29/15/. 


Uit x5e = Rs (cos. p + isin. ©) = Rz cos. @ + Rg sin. p X d 
volgt dan: 


COS. © = 


X5,6 = 2.8003 + 2.6454ï. 
Bij R, krijgt men als argumentvergelijkingen: 
cos.8 o — 1.69792 cos.* @ +- 0.92896 cos. © — 0.144123 SS 


en 
cos? 9 — 0,97998 cos. 9 + 0.186848 = 0. 


De G.G.D. der beide eerste leden, = 0 gesteld, geeft hier: 
0.03855 cos. © — 0.009979 = 0; 


9979 
COS. 9 = 38550 log. cos. © = 9418307 — 10; 
p= TANDS OG 


1) Met het oog op latere vermenigvuldigingen voor de berekening van Rê, Rö enz. 
zijn hier met kleinere cijfers nog 2 decimalen na de 5 eerste genoteerd. 


45 


%34 = Rj (COS. p + isin. @) = R‚ cos. 9 + Ry sin. Xi geeft dan: 
N34 = 0.5993 + 2.2362 DG 
De argumentvergelijkingen bij R‚ worden: 
COS.* 9 — 1.48463 cos? + 0.60929 cos. o — 0.042065 —= O 
en cos? © — 1.02907 cos. @ + 0.229037 = 0. 
De GGD. der beide eerste leden is hier: 
— 0.08855 cos. @ + 0.062276. 
Gelijk nul gesteld, levert zij: 
B, log. cos. 9 = 9.84709 — 10: 
o= 45°18/50”. 
X12 = Ry; (cos. p+isin. ©) = R,; cos. pt R‚ sin.o Xi wordt nu: 
Xi, = 1.394 + 1.4149 X ú 
Contrôle. De verg. werd gevormd uit de wortels: 
Kane 20 (Vlo HBr 040T7A 0 t; 
Oe OON AM 0 == 0 Ort 2007 
Kia= 14 tij2=14 + 1.41421 X ú. 5) 
Hierbij zouden behooren: 
log. R; =log. V(7.84+7) =log. V 14.84 = 0.585715 (vroeger ge- 
vonden: 0.58571 56). 
log. R‚ =log. 5.36 =0.36458; (vroeger gevonden: 0.36457 50). 
log. R‚ =log./3.96 =0.29885; (vroeger gevonden: 0.29885 co). 
De waarden voor log. R gevonden waren dus tot ongeveer in 
de laatste decimaal juist; die voor de termen der complexe 
wortels verschillen reeds in de 4e decimaal van de werkelijke 
waarden. De oorzaak van deze afwijking is gelegen in de vrij 
uitvoerige berekening voor het argument, waarin ook hoogere 
machten van den berekenden modulus optreden en bovendien 
op het laatst door aftrekking nog een decimaal uitvalt. Boven- 
dien wordt deze onnauwkeurigheid hier niet meer te niet gedaan 
door een worteltrekking met hoogen exponent zooals bij de 
berekening van den modulus. Om dus het argument en dan de 
termen van de complexe wortels eveneens in 5 decimalen nauw- 
keurig te vinden zal men de berekening voor het argument in 
bijv. 7 decimalen moeten doen. 


COS. 0 = 


1) De gevonden benaderingen wijzen duidelijk op de waarschijnlijkheid van de zeer 
bijzondere reëele waarden 2.8, 0.6 en 1.4. Door onderzoek vindt men dan verder 
gemakkelijk de imaginaire termen iV/7, iV5 en iV2; ons onderzoek beoogt echter 
niet de nasporing van zulke bijzondere gevallen. 
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S 9. Andere bepaling van het argument der complexe wortels. 

Zij o, het argument van y‚ in de eindvergelijking. Denken wij, 
dat de wortels der achtereenvolgende getransformeerden niet 
tegengesteld genomen zijn, dan zal het argument van dien wortel 
der voorlaatste getransformeerde, die met y‚ overeenkomt, ge- 
weest zijn: 5 of 5 Hr. 

Uit de verdubbeling van beide waarden ontstaat nl. als 
argument (tusschen O en 27) o,. 

Een substitutie, met slechts zeer ruwe berekening, wijst terstond 
uit welke van beide waarden bij die voorlaatste getransformeerde 
behooren. Men kan op dezelfde wijze langs de achtereenvolgende 
getransformeerden teruggaan naar de oorspronkelijke vergelijking. 
De substituties kunnen geheel in reëele grootheden uitgevoerd 
worden; bij vergelijkingen met reëele coëfficienten met behulp 
van (11) of (12). 

Heeft men in de berekeningen de wortels der getransformeerden 
wel tegengesteld genomen, dan kan men gemakkelijk tot het 
hierboven gedachte geval terugkeeren door aan de coëfficienten 
der getransformeerden eerst om den anderen het tegengestelde 
teeken te geven. 

Men kan echter ook het tegengestelde nemen van de wortels 
als een vermeerdering of vermindering van het argument met 7 
opvatten. Men kan dus het geheele systeem der getransformeerden 
onveranderd laten, als men in de achtereenvolgende vergelijkingen 
voor het argument substitueert: | 


BEONE DN TD 

Worden ín het voeger behandelde voorbeeld Il de opeenvol- 
gende moduli (van de eindvergelijking afgerekend) R,, Rs, Ro, R, 
en R‚ genoemd en de bijbehoorende #’s voorgesteld door 04, 93, 
Oa, @, En 9, dan zal Rs = VR,, Ry = VR3 enz. Men vindt dan 
achtereenvolgens: 

0, = 93°3459”, 0; = 43°12/30'.5, 09 = 68°23'44".75, 
0, = 5548/76 en op = 117°54/3”.8. 

Carvallo, die zegt dat Encke „méconnaît l'idée féconde de 
Pinventeur,” omdat hij „n’emprunte au calcul que la connaissance 
du module” leidt ook uit ya door achtereenvolgende substituties 
het argument van z, af. Hij denkt hiertoe, — evenals voor de 


€ 


TRE 


afleiding van den coëfficientenregel, — de gegeven vergelijking 
en de getransformeerden gebracht in den vorm: 
o(z?) + zF(z°) = 0. 

Substitueert men hierin voor z? de waarde van — y, uit de 
eindvergelijking, dan vindt men uit: 
rde) 

402) 

(den modulus en) het argument van den aan y‚ beantwoorden 
den wortel der voorlaatste getransformeerde; op dezelfde wijze 
kan men opklimmen tot de oorspronkelijke vergelijking. 

Bij vergelijkingen van betrekkelijk lagen graad, waarbij een 
niet onbelangrijk aantal transformaties noodig was, levert echter 
Encke's methode eenvoudiger resultaten dan die door achter- 
eenvolgende substituties. 

Wanneer men de achtereenvolgende getransformeerden der 
boven behandelde 6e-machts-vergelijking onderzocht met de 
substitutie: 


ZZ == 





fi en Pk 
Pk ) 25 


heeft men voor x56 achtereenvolgens: 
Pe 76° 1720", 0, =128°840", o,=154°420”, mg = 167°2107, 
ps = 6°28'55”, w, = 93°14/27".5, or = 43°22'46”". 
Hierbij vindt men dan: 
Rs cos. © = 2.79994 en Ry; sin. 9 = 2.64588, 


dus X56 = 2.79994 + 2.64588 X i. 


Deze oplossing stemt nauwkeuriger overeen met 
X56 = 2.8 + iV7 —= 2.80000 + 2.64575 X í 
dan de vroeger gevondene. 

Een oordeel, aan dit onderzoek door substitutie verbonden, is 
de omstandigheid, dat de achtereenvolgende @'s uit de vorige 
worden bepaald door telkens door 2 te deelen. De fout in het 
argument van y‚ wordt dan ook telkens door 2 gedeeld. 

In ons voorbeeld [V kan men zich bij deze substituties zeer 
beperken door op te merken, dat men voor @, @4, 03, En @ 
nauwkeurig genoeg handelt door alleen op de drie eerste coëtifi- 
cienten te letten. Deze drie coëff. kunnen binnen de grenzen 
eener ruwe benadering gelden als de coëff. der vierkantsverge- 
lijking, die x, en x, bevat. 
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Omdat de middelste coëff. dan —= —2R cos. 9 moet zijn, blijken 
Qs, P4 en 93 in het 2de of 3de Quadrant te moeten liggen; ov, zal 
weer in het eerste quadrant gelegen zijn. Voor 9, en ox moet 
men de substitutie inderdaad uitvoeren. 


V. De methode bij transcendente vergelijkingen. 


S 10. In de omgeving van een regulier punt 2 = a kan een 
functie f(z) volgens het theorema van Taylor—Cauchy worden 
ontwikkeld naar opklimmende machten van z — a. 

Deze ontwikkeling heeft een beperkt gebied van geldigheid 
(holomorphie). Voor alle binnen dit gebied gelegen waarden van 
z moet de waarde van f(z) met elken gewenschten graad van 
nauwkeurigheid kunnen worden gevonden uit de machtreeks: 


ECE 

n=0 
wanneer men & slechts groot genoeg neemt. Omgekeerd, als 
men een binnen het geldigheidsgebied gelegen waarde van z 
gevonden heeft, waarvoor de machtreeks = O wordt, zal ook f(z) 
met elken verlangden graad van nauwkeurigheid = 0 gemaakt 
kunnen worden. De nulpunten (wortels) van f(z) kunnen dus met 
elke verlangde praecisie gevonden worden door oplossing van 
een hoogere-machts-vergelijking, — dus door de methode van 
de scheiding door machtsverheffing. Als « niet = 0 is, behoeft 
men slechts z — a door y te vervangen om te komen tot: 

ECrpr=0. 

n=0 

Voorbeeld V. 1 + zsin.z= 0. 

Laguerre heeft deze vergelijking besproken in de Comptes 
Rendus *). Hij komt door overweging van het feit, dat de ont- 
wikkeling in het eerste lid aanvangt met 1 +2? en dus een 
lacune vertoont tusschen twee termeh met gelijk teeken, tot het 
vermoeden, dat de vergelijking imaginaire (complexe) wortels 
zal hebben. Dit vermoeden wordt dan tot zekerheid door onder- 
zoek van het „geslacht” van 1 + z sin. z, dat = 1 blijkt te zijn. 

Onze methode levert met zeer weinig moeite de waarde van. 
de beide kleinste wortels, die zuiver imaginair blijken te zijn. 


I) XCIV, 636. 
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Na ontwikkeling gaat f(z) = 0 over in: 
Bt ZN 
bien 

Vervangen wij eerst z? door y; het transformatie-tableau 
worät dan: 


Omid + Bnn ani 


1 ; 1 1 
h k 6 0 Ne 0 
o1 i & KO Eion en: 
3 280 302400 5040? 


Ps 1 nl. 0.24461 nl. 6.09235—10 
Ze 1 nl. 0.48919. 


In elke kolom is de transformatie zoo ver mogelijk voortgezet 
alvorens met de volgende kolom te beginnen. | 

Na drie transformaties is de eindtoestand reeds ingetreden 
voor den tweeden coëfficient. Wij vinden: 


mod. zi = Terdden’ log. 0.48919 = 0.93202. 

Het argument van y°, (in de eindvergelijking) is 7; op de 
vroeger aangegeven wijze wordt licht gevonden, dat voor y,* 
y‚* en y, eveneens het argument = r zal zijn. 

Dat van zie in 1 + zsin.z=0 zal dus + bor moeten zijn. 
Zoo blijkt dan: 


zi = + 0.93202i. 


Contrôle. Bij substitutie dezer waarden in de gegeven verge- 
lijking wordt het eerste lid = 0.0000002. 

Voorbeeld VI: Z— esin. zZz — M= 0. 

Merken wij allereerst op, dat deze Keplersche vergelijking altijd 
een reëele oplossing zal hebben, gelegen tusschen M en M + e. 
Het eerste lid wordt voor deze waarden van z nl. respectievelijk 
gelijk aan: 
| | — esin. M en e —e sin. (M + e). 

De eerste uitdrukking is negatief, de tweede positief *). De 
ontwikkeling van het eerste lid kan voorgesteld worden door: 


EA cos. 92) — M= 0 
rn stra 120 CCS z) + M = 0. 
Wordt de term met z° verwaarloosd, dan zal er bij substitutie 


_1) Men mag nl. aannemen O0 < M < z. 
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der uit het resteerende deel der vergelijking voor z, berekende 
5 

waarde een fout ontstaan koogstens = En $ 

Zij nu M= 27°31’5”.23 of = 0.4802818 en 


log. e = 9.3897262 — 10 2); 





5 
eh is dan 0.0002 ongeveer. 
3 
M (lez F=0 


heeft dan bij deze waarden voor M en e het volgende trans- 
formatie-tableau: 


l —num.log.0.19626 — num. log. (8.93007 — 10) 
jl 0.39252 —(9.42736 — 10) — (7.86014— 10) 
1 0.82158 9.03081 —10 5.72028— 10 

| 1.64103 8.03903 —10 

| 3.28205. 


Ook hier is elke kolom zoover mogelijk voortgezet alvorens 
met de volgende te beginnen. 

Wij vinden nu: 

16 1 
Zin Vlrumolon 3 28200 weeen 
In de hoekmaat ís dit: 
Zi = 90° 43/1677, 

Contrôle. Volgens Oppokzer is z= 35°43/30”.50. Tot deze 
waarde zouden wij veel dichter genaderd zijn, als de term met 
z? nog in rekening genomen was. De fout bij substitutie zou 
DDS di. ongeveer 0.000005 bedragen. 

Voor een eerste benadering, waarop men licht een bekende 
differentiaal-benadering kan laten volgen, is een nauwkeurigheid 
tot op minder dan 1 minuut echter meer dan voldoende. 

In de omgeving van een m-voudige pool a kan een functie f(z) 
nog steeds worden ontwikkeld in een machtreeks, maar deze zal 
dan ook negatieve machten van z— a bevatten en dus voor- 
gesteld moeten worden door: 


n=k 
2 Cole — a)”. 
m 


n= 





dan hoogstens = 








1) Voorbeeld uit Oppolzer, Bahnbestimmung der Plan, und Kom. 1, 56 met een 
kleine wijziging: M is naar het 1e quadrant overgebracht. 
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Het is echter gemakkelijk in te zien, dat ook hier de macht- 
reeks in plaats van de functie genomen kan worden voor het 
onderzoek naar wortels binnen het gebied van convergentie 
(holomorphie), als slechts # groot genoeg genomen wordt. Door 
vermenigvuldiging met (Z — a)” en vervanging van z — a door y 
krijgt men weer een gewone hoogere-machts-vergelijking. 


S 11. In de voorafgaande bladzijden zullen nu de voordeelen 
duidelijk genoeg geworden zijn van een methode, die door een- 
voud en algemeenheid alle andere ver achter zich laat. 

Er blijven nog wel enkele vragen ter beantwoording over, met 
name die naar de beste wijze van voortzetting eener eenmaal 
gevonden benadering. Als men een wortel bijv. tot in 10 decimalen 
nauwkeurig moet kennen, zou het een zeer tijdroovend werk 
zijn, alle berekeningen in dat aantal decimalen te verrichten. De 
bespreking van de beste wijze om zulke verfijnd-nauwkeurige 
resultaten te bereiken wordt bespaard voor een volgend hoofdstuk. 
Als uitkomst van dit onderzoek zij hier echter reeds medegedeeld, 
dat mi. de eenvoudigstede eenvoudigste wijze om aan een eerste 
benadering zo van een wortel eener vergelijking f(z) = 0 een 
correctie van willekeurig vastgestelde nauwkeurigheid aan te 
brengen bestaat in de ontwikkeling van dezen wortel z in de 
omgeving van zo als (implicite) functie van f(z). 


HI. VOORTZETTING DER BENADERING. 
UITBREIDING VAN DE METHODE VAN NEWTON., 


De benaderingsmethode van Newton voor de oplossing eener 
hoogere-machts-vergelijking f(x) =O onderstelt, dat men een 
benadering X9 voor zekeren wortel heeft gevonden, zoodat de 
werkelijke waarde kan voorgesteld worden door x= Xx + A, 
waarin A een kleine correctie voorstelt, die nog aan X, moet 
worden aangebracht. 

Deze A is bepaald door de voorwaarde: 

JAA EO LS RR (1) 
welke door ontwikkeling van het eerste lid kan worden ver- 
vangen door: 


INNER 
Liz (n—]l)! 
Voor de rigoureuse berekening van A zou dus eigenlijk weer 
de oplossing eener vergelijking van den nen graad noodig zijn. 
De omstandigheid, dat A als eene grootheid van gering bedrag 
gedacht wordt, kan echter groote vereenvoudiging opleveren. 
Het is overbekend, dat Newton om een eerste correctie te 
vinden alle termen met hoogere machten van A dan de eerste 


verwaarloost, een benadering van A berekent uit: 


_ fo) 
fo)’ 
dan X+ Ao =X, als benaderde waarde voor x beschouwt en 
met x, dezelfde redeneering herhaalt, als eerst met x) gehouden is. 
Even bekend is het, dat ook de voorwaarden, waaronder deze 
methode werkelijk tot het beoogde resultaat voert reeds door 
haren genialen ontdekker zijn ín- het licht gesteld. Door den 
Hornerschen Algorithmus is men verder ín staat gesteld de voor 
deze methode benoodigde berekeningen zeer elegant en met 
groot gemak uit te voeren. 


EE EA 2 TA, ATS OEE 





VAN Wen 





op) 


/ 


Het is zeer begrijpelijk, dat men bij deze methode liever enkele 
malen meer het formuletje 
fx) 

EUT UR EERE dE) 
toegepast, dan meerdere termen in de ontwikkeling (2) in rekening 
te brengen. Reeds door den term met A* in acht te nemen zou 
de berekening in het algemeen zooveel meer werk vereischen, dat 
het voordeel eener grootere nauwkeurigheid van de eerstvolgende 
benadering daardoor in de schaduw gesteld werd. Het spreekt van 
zelf, dat het toevoegen van een term met A® en termen met hoogere 
machten deze correctie tot groote complicaties zou voeren en weer de 
oplossing eener nieuwe hoogere-machts-vergelijking zou vereischen. 

In deze eenvoudige overweging ligt naar allen schijn een ge- 
heel voldoende verklaring opgesloten van het feit, dat de z.g. 
Newton-Hornersche methode, werkende met het formuletje (9), 
stereotyp van de eene generatie op de andere is overgegaan. 

Er zijn echter gevallen denkbaar, waarin het wenschelijk en 
door vereenvoudigende opmerkingen ook gemakkelijk uitvoerbaar 
is om meer termen der ontwikkeling (2) tot hun recht te 
doen komen. ' 

Een eerste opmerking is de volgende. 


Lade Le ) geheel weg 





In plaats van, evenals alle verdere termen, A? 


te laten zal grootere nauwkeurigheid bereikt moeten worden door 
in den derden term van (2) voor A te substitueeren de eerst- 
gevonden benadering Ag en dus A te berekenen uit: 


fe) + A foe) + ALP 0 Baren 
ME Aad PAPE) 1e | ve El =0) 


In het eerste geval verwaarloosd men nl. een oan Nm: 





(%) 
Oe 
in het tweede slechts: 

(ASG INGE 


of, als A, dicht bij de werkelijke A ligt, zoodat A = A, +ò, 


1) Wij onderstellen voorloopig /\ van zulk een dimensie, dat van termen met /\3 
enz. kan worden afgezien; dit zal in het algemeen altijd kunnen geschieden, als /\ 


ongeveer 





l 
1000 van het bedrag van x is. 
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slechts een grootheid, die met groote benadering wordt voor- 

gesteld door: 
fo) %) 

2d A TERDK 

Uit (4) leiden wij af: 

NV = — fo) nk f(X) Ô f”(%o) (5) 
FG) Sd), 1.270) vn 
leder, die de Newton-Hornersche methode slechts enkele malen 
heeft toegepast, zal terstond inzien, dat A’ op eenvoudige wijze 
te berekenen is. De gewone correctieterm behoeft slechts ver- 
minderd te worden met zijn quadraat, vermenigvuldigd met het 
guotient van den derden en den tweeden term *). 

Denken wij nu volgens de methode van Newton of Gräffe een 
eerste benadering gevonden tot op een of. enkele duizendste 
deelen van haar bedrag nauwkeurig, dan ‘zal in den regel het 


fo) 


gewone formuletje A, = EE) niet, maar wel de nauwkeuriger 


en dus van veel geringer bedrag is. 








correctie (5) den gezochten wortel kunnen opleveren met de 
gewone precisie in 7 decimalen. 

Wij zullen hier voorloopig geen toepassing der correctie (5) 
geven en ook liever afzien van het beschouwen van verdere 
termen in het eerste lid van (2). Het zal nl, uit de volgende 
bladzijden blijken, dat (5) kan afgeleid worden als een bijzonder 
geval van een ‘zeer algemeene ontwikkeling der correctie, die 
gevonden wordt door implicite differentiatie, maar ook op elemen- 
taire wijze kan worden berekend. 


Elementaire omkeering eener functie. 


le Voorbeeld. Als gegeven is: 
KI ae ea en Sinn rens (5 
vraagt men omgekeerd x te „ontwikkelen” in y: 
nat By ppt tytn) 
Oplossing. Door 1 met zich zelf te vermenigvuldigen vindt 
men zeer gemakkelijk: 
ye = KI 2E ZAP A AAE Tr OKER en 


1) 2e term beteekent hier: laatste getal van de 2e kolom (van rechts) enz. 
2) Aan y= 0 beantwoordt blijkens (1) altijd een waarde x = 0; men kan 
dus « wel weglaten. 
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Door de vermenigvuldiging van (2) met (1) uit te voeren blijkt dan: 
Oe Or Zer DCA (Aat Ir 2e Ib, oe 0 (8) 
Men kan steeds hoogere machten van y op deze wijze in x 
uitdrukken. Tusschen (1), (2) en (3) kan men x? en x? elimineeren: 


vAn rin a aten ne RE 


—y= KBK... 
Aes air (aaide Lijten a 
+ 
yy tryi=x ETET 


Wanneer x een klein getal is, zal bij benadering 
Rn METEN 

De 4e en hoogere machten van x zullen dan zeer weinig 
gewicht in de schaal leggen. 

Men kan de vermenigvuldiging met (1) steeds voortzetten en 
steeds verdere machten van x elimineeren. Men vindt dan steeds 
verdere termen van de ontwikkeling voor x. De volgende term 
wordt — y* enz. Wij hebben dus: 


oh An cia zb And an ma REA 
Contrôle, Het ís licht in te zien, dat y = TZ (meetk. reeks!) 
Hieruit volgt: 


y—yX=X en dus EEP 


Lr 
Bij ontwikkeling (deeling) blijkt terstond: 
En mr det Do mal Ata Ad 
Ile Voorbeeld. 





ME GN + 5 7 Baene (SUS Tee ks Wenen (0) 
Nu sy lEtE)e 

= x%1 AE) Et 0% 

nd Nn or is Eet 

y= xl zi | tE. 
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Nu kunnen x? en x? geëlimineerd worden: 





x5 x5 dek Den factor van y? vindt men 
Ke nee 
J 3115! 7! uit ri 
ss ce Kdl 13x7 De factor van yö wordt ge- 
IE hi 2.3! rpm 4 vorderen ci 
5! 2,3! 
3 en 3X5 1 , samen te tellen. 
zE 
B SRAL 
ET keen rap 
Men ziet hier terstond in, dat ter eliminatie van x? weer ver- 
meerderd zal moeten worden met Se re 5 STAR 
Àls x klein genoeg is, zal dan zeer benaderd: 
5) 
5 Rete 
id Dn EH. A MEEL 


Contrôle. Gegeven was: y =sin.x; dus zal x= bg. sin.y moeten 
zijn. Nu is 


bg. sin. y = Vri ee Vree 


1. me 
VEER 
3 


mph He Re 
_N.B. De ontwikkeling (10) geldt alleen voor beperkte waarden 
van y zoowel als van x. Dit wordt later toegelicht. 
HIL. Algemeen Voorbeeld. 
v=Ark tr Aart Agt orn 
Wanneer gegeven was: | 
Vm Ag A te Arti Agee Mene 
zou men kunnen stellen: 
Ver A= SA Kr Ae Ar 


en X dus kunnen ontwikkelen naar opklimmende machten van 
y of y —À. Deel beide leden van Ae door ie 


Lekt d Bags. MODE 0 


Door (12) met zich zelf te OEE krijgt men: 


BAL Ars, Ao? + 2AA, 
Fa Ott Sgt, EE 
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Uit (12) X (13) vindt men: ° 























IV ve SAo a SAo + SAiAs , 
te vee it has NDR (1E) 
De eliminatie van x° en x° gaat nu aldus: 
ene Ul Age Dr oe 
gee in nen me he. 
ld Ke Ar 2As* A51 2A Ass, 
la, TETN RS CRE EBT EAO 
Ee 2A,* —AAs ; , OAoS—3A,AoAs , 
A? A el Kee AE A KLP NONE 
y zhe ef snee Ge ERM 
A sd A a TEM: 
Om ook xt te elimineeren blijkt uit bovenstaande sommatie, 
dat de factor van el zal zijn: 
1 
A ?A,-A9*-—2A AgA3t6A0°—3A, AoAs AA, —5A,AzAs0A0® 
eend REA zZz IN 3 e 
1 
Wij hebben dus de ontwikkeling: 
BREA Aa 
A, AA, Aj A, 
APA, — 5AA,Ag + Dep (15) 
A nae 


Bovenstaande formule kan worden toegepast voor een zeer 
nauwkeurige voortzetting der benadering van een wortel eener 
algebraïsche of transcendente vergelijking. 

De beide voorwaarden, waaronder (15) goede uitkomsten 


geeft, — nl. x en Je beide van kleine volstrekte waarde, — 


A, 
worden tegelijkertijd vervuld door aan te nemen, dat men een 
eenigszins benaderende aanvangswaarde heeft bepaald }). 
Zij bijv. ter oplossing gegeven: 
MEA ge kar Va vor. 20 DER Bn IPRA 
Zij verder x, een eerste benadering voor x, zoodat x= %) + A 
of % — X = A gesteld kan worden. 


1) Dit is juist de onderstelling, waarvan ook uitgegaan moet worden bij de gewone 
Newton’sche methode. 
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Verminder nu de wortels met x,. Dit kan met behulp van den 
Horner'schen Algorithmus geschieden, er komt dan: 
0: B TB Art De Ate ien ESAT 
B, B 


ee B, NN 
BEU A +5 Alk: Te 


Het resultaat van de Re van bovenstaande functie kan 
onmiddellijk worden opgeschreven, door in (15)-A te vervangen 
door B en y door — B. Dit geeft: 





PE 
BABE SBS RDR 
B, B, B, Rl el 
IE Dn 5, B, 8 ae Sp B en ZA (17) 


Als benaderingsformule kan dus gegeven worden: 
BiBi ee El Bie 


B BBR EB en 
B, BB, En 
—(5'— Sp ptp) pe 9 enz BEE 


Het is terstond duidelijk, dat de benaderingsformule van Newton: 





fx) 
bm ki pre Mekel ee 19 
btn KE 
een zeer bijzonder geval is van (17). | 
Immers 
0 = F(x) = f(X + A) = fl) + f'(%) A + 
LE At EA AEN (20) 
zoodat JO) DIEET Wi ZOE 
Tae B, 
N= Ko — B,’ 


Vergelijking van (16) met (20) geeft ook onmiddellijk de ge- 
legenheid om formule (18) in plaats van in de coëfficienten der 
hefleide gegeven vergelijking uit te drukken in f(X), f(X) en 
in de verdere afgeleiden van f(x). 

Eerst mogen echter enkele toepassingen van (18) hier plaats 
vinden. 


B 
1) Uit deze ontwikkeling blijkt als voorwaarde, dat de volstrekte waarde van En 
1 


klein moet zijn; dit is juist dezelfde voorwaarde als bij de benaderingsmethode van 
Newton, want B, = f(x) en B, = f(x). 
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Je Wen Kiezen EE de vergelijking: 
—_Txt1=0, 


die door Lagrange is beschouwd in zijn Traité de la résolution 
des équations numériques, Chap. IV. Wij denken ons als ge- 
vonden, benaderde waarde voor een der wortels x,= — 3,05. 
De volledige Horner'sche Algorithmus voor dit geval kan dan 
door onmiddellijke vermenigvuldiging met het getal — 3,05 in 
geheel (niet cijfer voor cijfer) aldus worden voorgesteld: 











Del 0 7 +7 
Oa 0IS0 25e 7022098 
—305 _ +23025 _ — 0,022625 = B, 


805 18,605 


— 6,10 + 20,9075 = B, 








— 3,05 
— 9,15 = B. 
5 Bonen, 02202010: 
Nu is — B,“ 20,9075 0,001082148 - 
2 
B, Bo Sd sÀ en — 0,00000 05125. 





__B,B,? 20,9075 °° B? 
Uit de ontwikkeling (18) blijkt gemakkelijk, (B; = 1, B, = 0), 
dat geen verdere termen noodig zullen zijn, als x in 9 decimalen 


nauwkeurig berekend moet worden. 
Men vindt dan: 


X= — 9,05 + 0,00108 2660 = — 3,04891 7340 .... 


Opmerking. Een groot voordeel van deze methode is, dat de 
termen — B —ö pijen logarithmisch berekend kunnen worden. 

Als in het behandelde voorbeeld een nog grooter aantal nauw- 
keurige decimalen verlangd werd, zou men de eerste correctie 
5 ook door rechtstreeksche deeling kunnen berekenen. 

1 
Ie Toepassing: 
x8 — 2x — 21 =0. 
‘Wij zullen deze vergelijking geheel volledig oplossen. 
Zonder eenig verder onderzoek passen wij de methode van 


Gräffe toe: 
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1 0 01 — 21 
oid 42 441 441 
9 1 882 157437 194481 
98 ‘1 nl 566563 nl. 1038817 © nl 1057776 
24 1 1121888 2077631) 21.15552 
95 1 929.41840 41.55262 42.31104 
9e 1 4483635 83.10524 84.62208 
PRA: 5 SIDE 002370: 1x SENDER 
loor 2259650 050705: AGNNENEN 
lop As ee — 0.70057: x/ = 5.0184. 


De teekens van de wortels volgen onmiddellijk uit x, +X, +-X3=0; 


het blijkt, d 


at: % = — 1.0561; % = — 3.9623; X3 = 5.0184. 


De voortzetting der benadering en tevens de herhaling der 
contrôle geschiedt nu het eenvoudigst door den Algorithmus van 


Horner op 
vermindert 


de oorspronkelijke vergelijking toe te passen; men 
met de gevonden benaderingen, maar kort deze in 


tot 2 decimalen. De berekeningen zijn dan nog zeer eenvoudig; 
wie ze bezit, kan dan bovendien de tafels van Crelle toepassen, 


waarin de 
geven. Men 


producten van alle getallen van 3 cijfers staan opge- 
heeft dan: 


a) Voortzetting der benadering van x,: 


jl 


Men vind 























0 — 21 — 21 
— 1.06 1.1236 21.068984 
— 1.06 — 19.8764 0.068984 = Bo. 
— 1.06 2.2472 
— 2,12 — 17.6292 = B. 
— 1.06 
— 3.18 = B,. 
t, logarithmisch: 
nk —= + 0.003 9130 54; 
1 
2 
En = — 0.000 0027 62; 
4 
2 
X=— 1.06 — ln Ti = — 1.056 0897 08. 


1) Wanneer alleen de kleinste wortel verlangd werd, zou de bewerking hier gestaakt 


kunnen worden. 
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b) Voortzetting der benadering van x,: 


























l 0 — 21 — 21 
— 3.96 + 15.6816 + 21.060864 
— 3.96 — 9.9184 + _0.060864 = B, 
— 3.96 + 31.3632 
— 7.92 26.0448 = B. 
— 3.96 
— 11.88 = B, 
rn 0.002 3368 96; 
B, 
2 
Ne = + 0.000 0024 91; 
1 
Dong 
Kg = — ENT — 3.962 3344 05. 
c) Voortzetting der benadering van x4: 
1 0 — 21 — 21 
5.02 25.2004 21.086008 
5.02 4.2004 + 0.086008 = B, 
_5.02 50.4008 
10.04 54.6012 = B, 
5.02 
BS Ba 
ole — 0.001 5532 41; 
B, 
BB 
Be? Ene 0.000 0006 65; 
En B, BoB, 
Xa = 5.02 — BEET 5.018 4460 94. 


De nauwkeurigheid der gevonden uitkomsten blijkt uit het 


feit, dat weer &, + x3 + X3 = 0.000 0000 00. 


Als resultaat der voorafgaande beschouwingen moge hier nu 


reeds opgegeven worden: 


De aanbevelingswaardigste volledige oplossings- -methode voor 
algebraïsche en transcendente vergelijkingen wordt gevormd door 


de combinatie van: 


a) de methode van Gräffe voor de bepaling eener aanvankelijke 


benadering; 


b) den algorithmus van Hörner voor de berekening der coëf- 


MOLE Be Br De eNZ.: 
c) de formule (18). 
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Omkeering van functies door implicite differentiatie, 
Door in (18) B, te vervangen door f(x), B, door f'(x), Bs 
() : 
Tue) . @) 
door gy en in het algemeen B, door GET gaat deze formule 
over in: 
(22 en 
Íx) LU RE, [21 Je Ie 1) 
RO) ANT SEWER EE AI TON 


of, eenigszins uitgewerkt: 








X= Xp 














G en fo zi Ee He En fa ft lik Vis hi 
rem ap pa jn je 
fv sd III AAN f* 
dr sn Sn 
of 
De fs fes ier gf ed f? 
ven ap pn sa | 
vi fe ine sE 
En 107 5 Lln RO: 


Formule (21) geeft de (aanvangstermen der) omkeering van de 
functie f(x) = 0 in de (beperkte) omgeving van de benaderde 
oplossing %o. | 

Deze formule, die geheel een differentiaal cachet heeft, kan 
met behulp der differentiaalrekening op gansch andere wijze 
worden afgeleid. 

Alvorens deze methode te behandelen mag echter de beteekenis 
van de omkeering eener functie wel wat nader worden toegelicht. 

Het le lid van de vergelijking: 


x3 + 2x? —X— 1 =0, of in het algemeen f(x) = 0 


is een éénwaardige functie van de grootheid x. Met elke wille- 
keurige waarde van x stemt slechts één enkele waarde van het 
le lid overeen. | 
Omgekeerd kan x beschouwd worden als een functie van het 
le lid der vergelijking, maar dan als een meerwaardige (hier 
driewaardige) functie. 
Wat dit beteekent blijkt eenvoudig uit het volgende onderzoek. 


1) Bij f, f’, f'‘ enz. moet men overal x, als argument denken. 
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Als aan het eerste lid niet de waarde 0, maar de waarde 1 
wordt toegekend, voldoen aan f(x)!)= 1 de drie wortels: | 
Ml, Aj EN Kz == — 2e 
Wij zullen nu de waarde van het le lid van 1 af langzaam 

laten afnemen en wel eerst van 1 tot 0.9. 
Er ontstaat dan de vergelijking: 
XS HX — KA —1=0.9 of x8 + 2x2 — x— 1.9 = 0. 
Hieraan voldoen dan: 
AEK er OI eN Mg 0 7 2 
Tengevolge van de afneming van f(x) met 0.1, neemt x, af 
met 0.02; x, neemt toe met 0.05 en «x; weer af met 0.03. 
Aan f(x) = 0.8 beantwoorden: 
eg == 0,90, en. Amr 206 3 9). 


Bovenstaande uitkomsten zijn met nog enkele volgende ín het 
onderstaande tafeltje vereenigd: 


f(x) X % X%3 
Ì l —Ì — 2 
0.9 0.98 — 0.95 — 2.08 
0.8 0.97 — 0.90 — 2.06 
0.6 0.93 — 0.81 — 2.12 
0.4 0.89 — 0.73 — 2.16. 


Hij blijkt duidelijk, dat aan een vloeiende verandering van f(x) 
drie vloeiend varieerende functies x,, X3 en Xx beantwoorden. 

Men kan het verloop dier functies x,, X en Xx; gemakkelijk 
graphisch voorstellen. Men neemt bijvoorbeeld” de waarden 
l, 0.9, 0.8 enz. van f(x) als abscissen en eerst de daarbij be- 
hoorende waarden van x, als ordinaten. Het verloop van x, 
wordt dan graphisch door een (kromme) lijn voorgesteld. 

Daarna wordt op dezelfde wijze x, en vervolgens x3 behandeld. 

De drie ontstane krommen worden dan 0 takken van de 
functie x = 9(y) genoemd. 

Elk dezer takken kan nu als (grafiek van) een gewone een- 
waardige functie beschouwd worden. 

Elk der wortels moet dus als een functie van f beschouwd en 
dan ook ontwikkeld kunnen worden. 


1) nl. aan x3 + 2x2 —x—l=l of x3 + 2x2 —x—2=0. 
2) Door de afkorting in 2 decimalen kan x, + X9 + X3 wel eens 0.01 van — 2 
afwijken. 
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Beschouwen wij nu eens de vergelijking: 
Je) =y. EE 
Fixeeren wij-in (22) een bepaalde waarde voor y en onder- 
stellen wij, dat x, een waarde voor x is, waarvoor f(x) de waarde 
Yo (dicht bij y gelegen) aanneemt. Nu zal x = v(y) elk der takken 
van de meerwaardige functie z (y) kunnen voorstellen. Wij kunnen 


























nu stellen: 
KK Ar en y= Vo tT Ayse terwijl f(X) Sn 
Dan zal 
BAO ee: EEE 
| dy. dyn FE 
dx 
He nel aA Ì Mik 1 Li 
| EN DR POET PN 
PT 
dx 
Ni d Ke MAS. d hit ‚| bie LFT TSP ET 
(y)= ls fi el ad dxf'? Afd nt f7 Gi Fi: ’ 
Iv E: ge Rr e, 1 hj SLOTS III en Er (3 112 mel) On 
(y) = sil Ti El f/n 
ffN— 1OF!f!f!!! + 15e 
gils Fi 
Uit So Piere (Mat A valt nu af te Res 
x= (Y) FP (Ho) Ay tr W)5 5 Le ) 5 nt ………… (24) 


Voor @(yo) kan % ERN Be vinden dan: 


Ji A? EN 3f’2 dr ERA A3) 


A 
x= Xt 5e 











HINA FS 
ffN LOF ff! + | 15f/3 Af, 
Ran fr 7 mat. Aan (25) 


Wanneer wij in plaats van f(x) =y de vergelijking f(x) = 0 
beschouwen, wordt blijkens (23) Ay = — Yo = — flXo) en (24) 

















wordt dan: 
is f rede is Oh foe A 
XS XT gk FLR gd 3E ak 
at En Hib) gl DCD 


Hier is formule (21) van blz. 62 door implicite differentiatie afgeleid. 


1) Men moet er steeds aan denken, dat f, f, f” enz. hier afkortingen zijn van 


fleeo), Fa), f'(x) enz. 
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Beknopte formule. 
De ontwikkelings-formule voor f(x) = 0 kan aldus worden ge- 
schreven: (zie 19): 


orbi 21 




















ad ai | 
Hierin is dan de iten De GE en ne 
CDO CDC XO 
D.C 
Mene aleemeen is.Cr = DCi XC = ae 


Mumkanook Cr DGX C= ee gesteld worden, als C‚ = x, 


genomen wordt. 
Het is dan licht in te zien, dat (21) aldus in beknopten vorm 
kan worden voorgesteld: 





id k 
Jo)* DoC, Co = X. Per (22) 


OE Of 

Formule (22) stelt in beknopten vorm de geheele, tot in het 
oneindige voortgaande, omkeering van f(x) = 0 voor in de om- 
geving eener benaderde oplossing x. 

Opmerking. Alle termen der ontwikkeling (zie 21) bevatten in 
den noemer f’. Hieruit kan worden afgeleid, dat ter berekening 
van x de aanvangswaarde Xx, zoo gekozen moet worden, dat 
deze dichter bij x gelegen is dan bij den dichtst nabij liggenden 
wortel van f'(x) = 0. 

Transcendente vergelijkingen. 

De formule (21) is zeer geschikt voor de oplossing door be- 
nadering van transcendente vergelijkingen. 

le Voorbeeld. x—sin.x=l of X— sin.x— 1 = 0. 

Men vindt gemakkelijk als benaderde KR LOSLDE, 

Bn kklen > IrS Loren) 
sin X — 0.9335 804 


Xo — Sin-X) = 1.0037 351 en dus f, = 0.0037 351. 





1) Men moet er steeds aan denken, dat f, f’, f'’ enz. hier afkortingen zijn van 


Á(xo), (0) Eee) enz. 
d 
2) D is geen factor, maar een operatie-teeken; DC, beteekent zeer: 


3) Zie bijv. Dr. B. Gonggrijp, Log. en gon. tafels en bijtafels blz. 145 en 146. 
5 


66 


Nu is f/= 1 —cos.x=1.3583 679; f/= sin. x= 0.9335 804; 
f”'= cos.x0 = — 0.3583 679; fIV = — sin.xo enz. 
Bij logarithmische berekening in 7 decimalen vindt men: 


L — 0.0027 49697; { eb ge 0.0000 0259 8208. 
iÈ NRJ 

Het is gemakkelijk in te zien, dat voor nauwkeurigheid in 
7 decimalen geen verdere termen in het 2e lid van (21) ver- 
eischt worden. 

Men vindt dan: 

X=% —/ — zi: —= 1.934 5632. 

Wanneer men deze waarde weer in graden, minuten en sec. 
uitdrukt, vindt men x= 1.934 5632 — 110°50’32”.25 1). 

Contrôle. sin: 140-5032725: == 934 SGL 

Ile Voorbeeld. Ke al 0: 

Omdat 3% = 27 en 4* = 256, zal er een positieve wortel gelegen 
zijn tusschen 3 en 4, 

Men kan gemakkelijk nog een decimaal bepalen. Beschouwen 











wij bijv. KON Kin ZR REE NLO AN 
Nu is 3.5 log. 3.9 = 3.9 X 0.544....=1:904.., 5 
9.6 108:3,6— 3:6,-0090r sat ON: 


Xo=d.0 is te klein, 09.6 te sroot; er zal echter eén waarss 
van x dicht bij 3.6 gelegen zijn. ; 
Wij kiezen dus x,—=3.6 en zullen den „Neperiaanschen log. 
van beide leden der gegeven vergelijking nemen: 
fo) =xlg.x—lg. 100=0 of 
fx) =xlg.x—4.60517 01859 88091. 


Niuzal fx) PD Oe f")=5 enz. 


Verder is Ig. 3.6 = 1.28093 38454 62064; hieruit volgt: 
X,lg.X0=3.61g. 3.6 = 4.61136 18436 63430 
lg. 100 = 4.60517 01859 88019 


fo = 0.00619 16576 75411; 
f’ = 2.28093 38454 62064; 


gene se ids 
Pie . TS Be 


1) Zie bijv. Dr. B. Gonggrijp, Log. en gon. tafels en bijtafels, blz. 145 en 146. 
2) Zie bijv. Dr. B. Gonggrijp, Log. en gon. tafels en bijtafels en hiervan tafel IV. 
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Men ziet terstond in, dat de omkeerings-formule niet verder 
behoeft te worden voortgezet dan tot: 
l 
acf 
RRA Or Oren Aerden 
B ZON 


Logarithmisch (Brigg. in 7 dec.) blijkt 


X == Xo 


El — 0.00271 4528 


Aal 
72 0.00000 0448 68 .... 
+ 


—=0.0027L 4977. 
Dus zal x= x, — 0.00271 4977 = 3.59728 5023. 


Opmerking. Deze uitkomst moet tot op ek nauwkeurig zijn. 


In Crelle XII is deze vergelijking op geheel andere en veel uit- 
voerigere wijze behandeld door M. A. Stern, (blz. 25 sqq.). Hij 
geeft als uitkomst op: x= 3.09728 5; hij vermeld verder, dat 
Maernmheelt afgeleid. — 3.59728 52 in” de Inst.” Calc diff., 
Bd. XXII, 1.—. 

In plaats van door „probeeren” nog een eerste decimaal te 
bepalen nadat het aantal geheelen van Xo gevonden is, kan men 
ook van de allereerste, ruwe, benadering x, = 3 of = 4 uitgaan 
en een eerste rekening met formule (21) in een beperkt aantal 
(bijv. 5) decimalen uitvoeren. 


Nemen wij bijv. x= 4; dan is Ig. x = 1.38629,; 


A0 054518 — 460517 —= 0,94001; — f/= 238629; Ff’ =— 
| PELL DE l Ur 1 
Nu wordt: 5 — 0.39393; 


1 fQ 
5e —= 0.00813; de volgende termen worden: 
(zie form. 21). 
……. = 0.00060 


………… = 0.00006 


BA 
A 





ale 
0.40232. 


68 


Dus vindt men als benadering: 
x= 4 — 040232 = 3.59728. 


Deze waarde is reeds tot in de laatste decimaal juist. Als men 
haar als uitgang neemt voor een nieuwe toepassing van (21), 
blijkt het, dat men dan reeds volstaan kan met 
B 

Fr 
Men vindt dan de boven reeds opgegeven waarde 

Xx = 3.99728 5023 terug. 


X == Xp 


IV. DE KEPLER'SCHE VERGELIJKING. 
E—esin.E = M. 


Reeds meer dan 300 jaar geleden stelde, Kepler het schijnbaar 
eenvoudige probleem op: 

Uit een gegeven punt van de middellijn van een gegeven cirkel 
wenscht men naar den omtrek een rechte te trekken, die het 
oppervlak van den halven cirkel, waarin zij gelegen is, verdeelt 
in de verhouding mm: n. 

Hij was er in geslaagd de volledige oplossing der elliptische 
beweging der planeten 
om de zon tot boven- 
genoemd vraagstuk als 
kern te reduceeren. 

Deze planimetrische 
kern is een zuivere twee- 
lingzuster van de Kep- 
ler'sche (transcendente) 
vergelijking. 

De laatste. kan op de 
volgende wijze uit de 
eerste worden afgeleid. 
(Zie de fig.) 

Zij F een gegeven punt van de middellijn AP van den gegeven 
(halven) cirkel AGP. Nu wordt verlangd FG zoo te trekken, dat 
OppsPRGROppan AFG! np: n 

Roem OMP =E, Nuis 

1 


Opp. PFG = Opp. PMG — Opp. FMG = ME, RSE 
=}R?.,E-—4R.MF sin. E. 





Verder zal 
Opp. AFG=Opp. AMG Opp.FMG= Es Es +5 MEER Sin, Be 
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Verlangd wordt dus: 


REE RMP sin. B): (TR R.MFsin.E)=min of 


5 Rr (mtn)=(t RER. MF sin. B): m, dus 
mrR? =(m + n)R?.E — (mm + n)R. MF sin. E. 
Na deeling van beide leden door (me + 7) R° komt er: 


Bie sin Be 
E — esin. Ei M. 

Het issdiudelijk;sdats0 SS Zinnen s OA Nen 

Het planimetrische probleem en de transcendente vergelijking 
van Kepler komen beide neer op de bepaling van de „exentrische 
anomalie” E uit de „middelbare anomalie” M bij gegeven excen- 
triciteit e. 

Meetkundig en algebraïsch hebben vele groote wiskundigen 
hunne krachten aan Kepler's probleem beproefd. 

Zuiver graphisch is het probleem niet (alleen met passer en 
lineaal) op te lossen. Wel gelukt dit, wanneer bepaalde kromme 
lijnen, zooals een sinusoide of cycloiden als zuiver te construeeren 
krommen worden toegelaten. 

Reeds in 1659 schreef J. Wallis in zijn werk „De Cycloide” 
zijn artikel „de problemate Kepleriano per cycloidem solvendo.” 
Hierin wordt het geometrische probleem, zooals het oorspronkelijk 
door Kepler gesteld was opgelost door middel eener „verkorte” &) 
cycloïde. De oplossing en het bewijs van de juistheid der con- 
structie is zeer omslachtig. Algebraïsche toepassing voor de nauw- 
keurige oplossing van Kepler's vergelijking blijven achterwege. 

Het ware nut van de graphische oplossing van het meetkundige 
probleem ligt echter in het eenvoudig en spoedig verschaffen van 
een benaderde oplossing van de Kepler'sche vergelijking, die daarna 
zoo nauwkeurig als men zelf verkiest kan ontwikkeld worden. 

Spoed en gemak zijn de twee omstandigheden, waardoor zulk 
een graphisch-analytische methode de voorkeur verdient boven 


of 





1) De namen „verkorte” en „verlengde’’ cycloiden zijn dikwijls omgewisseld. Bij 
Wallis heet de door mij „verkort” genoemde cycloide juist de „protracta.’’ Het komt 
mij echter voor, dat een cycloide, welker totale lengte kleiner is dan die der „zuivere 
cycloide verkort genoemd moet worden. 
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de vele reeds voorgestelde oplossingen, ook boven de beroemde 
analytische methode van LAGRANGE. 

Deze laatste moge hier eerst even toegelicht worden. 

LAGRANGE’s omkeeringsformule geldt voor een algebraïsche of 
transcendente vergelijking, die in den vorm gebracht is: 

X — ap (£) = Ko Of 
X— Xo 
RTD 

waarin a een geringe waarde heeft en x, een benaderde oplossing 
voorstelt. 

In een beperkt gebied van convergentie geldt dan de ont- 
wikkelingsreeks: 





oeren aes late ate bf (1) 





Boe ge 
x= Kot E De, 19 0) $7 Rien (2) 
De vergelijking van Kepler: E —esin.E =M of 
E —M 
sin. Po RK =d 0 aaneen e (3) 


is als het ware pasklaar voor de toepassing van de reeks van 
LAGRANGE. 
Deze geeft de zeer elegante uitkomst: 


jen 
E-M+È zr Dr (BEREN ek ed (4) 


Hiermede schijnt Kepler's vergelijking voor goed en onver- 
beterlijk opgelost. Maar, — om deze oplossing voor becijfering 
„aan te wenden moet men sin.”*M expliciet in sin. M. en in de 
sinussen van de veelvouden van M uitdrukken en daarna de 
differentiaties uitvoeren. Dan vindt men: 


2 8 
EM +esin.M +“ sin. 2M+(3 sin. 3M — sin. M) + 


4 5 - 
+ s sin. AM —sin.2M)+aaz (125sin.M-—81sin.3M+2sin.M)}-…. (5) 
Over de voorwaarde van convergentie dezer reeks 
Er NOEZ CAL, ) 


schreven o.a. LAPLACE en later STIELTJES. 
De elegante formule (4) blijkt in de praktijk grooten omslag 





1) In mijn verhandeling over nulpunten en oneindigheidspunten van 1908 komen 
verdere beschouwingen voor over de reeksen van LAPLACE en BüRMANN. 
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te vereischen, vooral als e een vrij belangrijke waarde heeft, bijv. 
e=0.5 ...., hetgeen bij de banen der z.g. kleine planeten 
dikwijls het geval is. 

Om deze reden maakt men in de Astronomie dan ook weinig 
gebruik van de reeks van LAGRANGE; men zoekt liever òf door. 
graphische middelen of door de methode der achtereenvolgende 
benaderingen een vrij nauwkeurige aanvangswaarde en brengt 
dan de noodige correctie aan door de Newton'sche benadering 
een- of enkele malen toe te passen. 

De methode, die hier zal worden gegeven, gaat ook uit van 
een graphische of numerische aanvangswaarde, maar dan gevolgd 
door de in het vorige hoofdstuk behandelde impliciete ontwik- 
keling (of uitbreiding van de methode van Newton). 


A. Graphische Benaderingen. 
L De verkorte Cycloïde. 


In 1901 verscheen fmijne graphische methode met behulp der 
verkorte Cycloïde in de Astronomische Nachrichten (Bd. 155, 
blz. 369). Eerst na deze methode ontwikkeld te hebben kwam 
Wallis’ opstel te mijner kennis. Trouwens zijn meetkundige be- 
schouwing is geheel iets anders dan mijn graphische oplossing 
der transcendente vergelijking. Deze is de volgende. 





Men denke zich een cirkel met straal M‚O =r, die over de 
rechte basis OV rolt. Het punt O zal een zuivere cycloïde be- 
schrijven; een punt P, op den afstand MoP van het middelpunt M,‚ 
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gelegen, zal een verkorte cycloïde beschrijven, aangegeven door 
DD 

Als nu O, een willekeurig punt van de zuivere cycloïde is, 
wanneer verder M, de plaats van het aan O, beantwoordende 
middelpunt ©) en < O,M,F = E de hoek van wenteling is 
(M‚F L OV), dan zal de abscis OF == rE zijn. 

Het punt P, dat de verkorte cycloïde beschrijft, zal gekomen 
zijn in P,, op de lijn M‚O, zóó dat MP, = MP = b. 

De abscis van P, zal zijn: 

x= OF — KF = rE — b sin. E. BIDEN le). 
Vergelijken wij met (6) nu: | 
M=E—esin. E of Mr = rE — er sin. E, 
dan blijkt er volkomen overeenstemming tusschen beide verge- 
lijkingen te bestaan, wanneer 
| ber en x= Mr 
genomen wordt. 

Hieruit volgt als grapkische oplossing der Kepler’'sche ver- 
gelijking: 

Men construeert de verkorte cycloïde, die beschreven wordt 
door een punt, dat gelegen is op den afstand er van het middel- 
punt van een over een rechte lijn rollenden cirkel met straal 7. Op 
de abscis zet men de lengte Mr = OK af en bepaalt het snijpunt 
P, van de ordinaat van K met de kromme. Uit P, als middel- 
punt beschrijft men een cirkel met PM, = er als straal, die de 
rechte M‚M snijdt in M, en laat de ordinaat M‚F op de abscis 
neer. Dan is < P,M,F de gevraagde hoek E. 

Opmerking 1. Men kan den straal van den rollenden cirkel 
‚(de generatrix) = 1 stellen; OK wordt dan = M. Wanneer OV en 
M‚M met- een verdeeling in radialen of graden (of in beide!) 
voorzien worden, kan E=OF=M M, afgelezen worden. Wanneer 
men een teekening op millimeter-papier heeft, kunnen beide 
ordinaten KP, en FM, achterwege blijven. 

Opmerking IL. Zelfs de cirkel met P, als middelpunt en e als 
straal kan achterwege gelaten worden. immers P,D = ecos. E 
en KP, = 1 — e cos. E; hieruit kan men E gemakkelijk bepalen. 
De aflezing van 1 —ecos.E=KP, is nog van groot belang 


1) Dit middelpunt beschrijft de rechte M‚M. 
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omdat men daardoor behalve de excentrische anomalie ook den 
voerstraal p van de planeet vindt. Deze is nl. 
p=a(l —ecos.E), (a is de halve groote as), 
zoodat men heeft: 
£ — y of En KE 

Ook de ware anomalie v kan door een eenvoudige constructie 
bepaald worden (zie fig. M‚,O,BA of MRSW). 

Opmerking MIL. De gegeven graphische oplossing schijnt aan 
een groot bezwaar onderhevig te zijn, wanneer men tot numerische 
oplossing van de Kepler'sche vergelijking wenscht over te gaan. 
Daartoe moet namelijk een verkorte cycloïde vrij nauwkeurig 
geconstrueerd kunnen worden. 

In het genoemde artikel van de Astronomische Nachrichten 
heb ik echter aangetoond, dat zulke krommen op mechanische 
wijze en zeer nauwkeurig kunnen worden geconstrueerd. 

Ik heb verder aangegeven hoe een graphisch tableau voor 
het geheele planeten-systeem (dus voor veranderlijke waarden 
van €) kan worden ingericht, zoodat men door één blik op dat 
tableau de aan zekere gegeven waarde van M beantwoordende 
waarde van E‚ — en ook de andere grootheden, — kan aflezen. 

Deze uiteenzettingen zouden echter te ver vallen buiten het 
kader der hoogere algebra. Wie er belang ín stelt kan ze ín 
genoemd artikel nagaan; misschien gaat een practicus eerstdaags 
het aangegeven idee verwezenlijken en daarmede een dienst be- 
wijzen aan alle sterrekundige observatoria. 

Wanneer een enkele aflezing E reeds in 3 min. &) nauwkeurig 
opleverde, zou een enkele Newton'sche correctie of liever de _ 
vroeger beschreven uitbreiding daarvan de grootheid terstond in 
honderdsten van. secunden opleveren. 

Zoolang echter het „tableau’” nog geen werkelijkheid geworden 
is, kunnen beide volgende methoden dienst doen. 


U. De Sinusoïde. 


In de Astronomische Nachrichten, Bd. 69, blz. 177, heeft 
prof. Dubois de volgende graphische methode voorgesteld. 


1) De implicite ontwikkeling kan reeds bij veel geringere nauwkeurigheid der aan- 
vangswaarde het gewenschte resultaat leveren. 
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Als men in de vergelijking: 


X—esin.x=M 
stelt: 
MENS, a Aat a Pe ek 
heeft men daarbij: 
X—ey=M Rt Annee (0) 


Nu is (7) de vergelijking eener sinusoïde en (8) die eener 
rechte lijn. Om x te vinden behoeft men slechts het snijpunt van 
(7) en (8) te bepalen en daarvan de abscis af te lezen. 

Men zie Aanhangsel 1. 

Prof. Dubois stelde voor den hellingshoek van (8) en daarna 
de lijn uit het punt x= M (y= 0) met dien geconstrueerden 
hellingshoek te trekken. Eenvoudiger komt het mij echter voor 
de lijn (8) te construeeren door de beide volgende punten te 
verbinden: 


Het eerste punt ligt op de X-as (Om), het tweede op de 
horizontale lijn y = 1 (Mr). 

In de figuur van Aanhangsel I moet M evenals e in absolute 
maat worden uitgedrukt. 


le Voorbeeld: 
0204 sin Ei 27 1314472, 


Nu is M= 0.475 0918 en M + e= 0.842 6142. 

Korten wij M en M +e tot in 3 decimalen in, dan wordt in 
de figuur OA = 47.5 mM. en MB = 84.3 mM. genomen. De lijn 
AB snijdt de sinusoïde in P. Als abscis van’ P lezen wij af 
71.6 mM. en kiezen daarom x= 0.716 = 41°. 


Opmerking 1. De sinusoïde ís hier geheel in absolute maat, 
d.w.z. ordinaten en abscissen zijn in den zelfden maatstaf ge- 


teekend. Bijv. bij x = 30° = 2 = 0.52360 (= 52.36 mM.) is 


y= sin. X = 0.50000 (= 50 mM.). 

Men kan ook elke mM. der abscissen-as een graad laten 
voorstellen; 10° wordt dan aangegeven door een abscis = 10 mM. 
enz. Dan wordt verder y= sin. 10° =0.17453 wel door een ordinaat 
—= 17.453 mM. aangegeven. De ordinaten zijn dan alle gelijk aan 
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die in Aanhangsel 1, maar de abscissen zijn ingekort, nl. ver- 
180 
31409 
Men moet dan om de lijn te construeeren, die de punten 
x=M ndi 2 
ans epi | 
verbindt ®), e ín graden enz. uitdrukken. Men vindt x dan onmid- 
dellijk ook in graden enz. 
2e Voorbeeld. Zie Aanhangsel II. 
E — 0.197 4769 sin. E = 131°41/18”.07, 


menigvuldigd met 


Oplossing. Bh dl 
AA OL 
Mate =ládn 


Men heeft dus de twee punten te verbinden: 
ah BNA Ol dS 8 
Hm | el 
Koren Ee 
Opmerking Il. Men kan desverkiezende bij het aanbrengen van 
de lijn x— esin.x= M altijd links van M= 90° (dus in het 
eerste quadrant) blijven. 
Is nl. M > 90°, dan kan men stellen: 
x= 180° — x/ en M= 180° — M’, 
„De verg. gaat dan over in: 
180° — x/ — e sin. (180° — x/) = 180° — M/ of 
Ke ae SIR AN 
Bij y = sin. x/ heeft men de rechte lijn bepaald door de punten: 
Kuis INK N= M’ —e 
y =0 ki y=l 
Men vindt dan £/ en behoeft nog slechts x= 180° — x 
te nemen. 
3e Voorbeeld. Zie Aanhangsel 1. 
E — 0.413 2767 sin. E = 126°47/7”.18. 


Mis 2120212 N00 
e = 0.413 2767 e = 0.41328 


M +e = 2.626 0979 MW — e = 0.51549 


| en vindt dan: 


/ 











1) De lijn x— ey = M behoeft niet geconstrueerd te worden; men kan gewoon 
een lineaal langs de beide punten leggen en leest de abscis alleen af van het snijpunt 
van de lineaal met de sinusoïde. 


AT 


Op de eerste manier krijgt men E=2.47 .... (Zie de lijn CQD); 
op de tweede vindt men E' == 0.67.,.. (Zie de lijn ERF). 
NMS delier en 0007 oe ZA, 


IL. Drie systemen rechte lijnen *). 


Stel in E —esin.E=M 
Hire en VM; 
er ontstaat dan: E —xsin.E =y. 


Op deze wijze ontstaan drie systemen van rechte lijnen. 
In Aanhangsel III zijn de ordinaten in graden uitgedrukt en 
elke graad is voor de eenvoudigheid = 1 mM. genomen. Een 


abstis als x — 0.193... (of e= 0.193...) is voorgesteld door 
19.3... mM. enz. De lijnen E —xsin, E=y zijn als volgt 
geconstrueerd: 


Als bijzondere punten zijn genomen: 
Kri 0 Onl 
AE F WEES 
De waarde E — sin. E‚ op de lijn x= 1 afte zetten, is telkens 
weer in graden uitgedrukt om den maatstaf op deze lijn x= 1 


precies even groot te laten zijn als op x = 0 (Y-as). 
4e Voorbeeld. 
E — 0.245 3162 sin. E = 27°31’5”.23. 

Zie Aanhangsel III. 

Men zoekt het snijpunt van x= 245 mM. (// Y-as) en 
y= 27°.5 (// X-as). Dit blijkt te liggen tusschen de schuine lijnen, 
die bij E= 35° en bij E = 40° behooren. Op het oog valt nu 
ligt te schatten: E‚= 36°. 

Aanhangsel III geeft een zeer beknopt algemeen graphisch 
„tableau”’ van de oplossing der Kepler’sche vergelijking. Het laat 
echter geen scherpe interpolatie toe, zoodat het mij nog steeds 
wenschelijk voorkomt, dat het cycloïden-tableau tot stand ge- 
bracht wordt. 

Opmerking. Het is hier wellicht de plaats om te laten zien, 
dat men geen waarden van M > in toepassing behoeft te 
brengen. | 


1) Deze methode is het eerst besproken door Rapau (Bulletin Astr. 1 382 en 
XI 289). 
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Als nl. M > zr was, zou men kunnen stellen: 
M = 2 —M’ en E = 2 —E’. 
Dan zou 360° — E/ + esin. E/ = 360° — M’ 
of E’ — esin. E’/ = M. 
De vergelijking is dan onmiddellijk herleid tot een andere van 
geheel denzelfden vorm, waarin het 2e lid < 7 is. 


B. Numerische benaderingen. 


L Uit x=esin.x— M blijkt, dat bij kleine waarden van e, 
de grootheid x weinig van M zal verschillen. 
Als eerste benadering kan men dan stellen: 
Ne ve oi. ANN 
hierop kan men als verdere benadering laten volgen: 
Ko ERSTE + M; 
vervolgens Ks CST NET 
Door LAURENT is bewezen !), dat de waarden %,, .%, enz. 
werkelijk naar de waarde van x convergeeren, als e <1. 
ge Voorbeeld. (Zie 4e Voorbeeld.) | 
x— 0.245 3162 sin. x — 27°31’5”.28. 
Als, ter vereenvoudiging, M = 27°31’ gesteld wordt, vindt men 
Ki 341 S= IO ZO ND 
IL. Voor een eerste benadering kan men bij betrekkelijk lage 
waarden van e ook wel 3 of 4 termen van de ontwikkelings- 
reeks van LAGRANGE nemen. 
6e Voorbeeld. (Zie 1e Voorbeeld.) 
X — 0.367 5224 sin. x = 27°13/14”.7. 
Fer TiS: 
M = 0.47509 
e sin. M = 0.16811 


2 
Dr sin. 2M = 0.05495 


Xo = 0.69915 —= 40°4’. 
N.B. Wij vonden vroeger graphisch met de sinusoïde 
Kir Deken 


1) Dit bewijs is o.a. te vinden in mijn verhandeling over „Nulpunten en Oneindig- 
heidspunten”’ van 1908. N 
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De benadering, welke 3 termen van LAGRANGE's ontwikkeling 
geven is van beteekenis minder nauwkeurig. De volgende term: 


eije 
en (3 sin. 3M — sin. M) 

geeft nog: 0.01558. Dit bedrag, gevoegd bij 0.69915 levert op: 
Xj = 0.71473 = 40°57’. 


Bovengenoemde feiten zijn een sterke aanbeveling voor de 
graphische benaderings-methode. 


UL De Regula Falsi. 


Het eerste lid van: 
Xx —esin.x—M=0 

wordt voor x= M. gelijk aan — esin. M en 

5 MANO „ e— esin. (M + e). 

De eerste waarde is zeker negatief, omdat OS MS r en de 
tweede is beslist positief. 

Er moet dus een wortel zijn tusschen M en M +e. 

Als men de werkelijke waarde van den wortel M + d noemt, 
stelt men volgens de regula falsi de evenredigheid op: 

dre=esin M: je —e sin. (M + e) +e sin. M$. 

(De toename van de functie wordt evenredig gesteld aan de 
toename van het argument.) 

Hieruit volgt dan: 

ik esin. M 
“T1 sin. M — sin. (M + e) 





e sin. M 
Aes LRE 1E sins Moes sinm0Mbe he) At MEN 9) 


esin. M 
— 2 cos. (M + de) sin. Je’ 


Bij kleine waarden van e zou ook kunnen: 
esin. M 


SMF Te cos. MF 49) 
6e Voorbeeld. (Zie 3e Voorbeeld) 
X — 0.413 2767 sin. x = 126°47'7.13. 
Berekening volgens (9) in 5 decimalen geeft: 
X = 2.46589. 
Esch vonden wij vroeger: E‚= 2.47. 





of X=M+t3 
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De, — voor zoover mij bekend is, — tot nu toe nog nimmer 
vertoonde en toegepaste eenvoudige formule (9) of (10) verdient 
als benaderingsformule verre de voorkeur boven de beroemde 
ontwikkeling van LAGRANGE. 


IV. De Methode van Gräffe. 


De algemeene oplossings-methode eener transcendente verge- 
lijking werd op blz. 48 ín S 10 besproken. 

Het eerste lid wordt ontwikkeld naar opklimmende machten 
van het argument. Zonder eenig voorafgaand onderzoek wordt 
de bekende transformatie toegepast. De bewerking wordt met 
elke kolom zoover mogelijk voortgezet alvorens de rekening met 
een volgende te beginnen. De overgang van een term van zekere 
kolom in zijn quadraat geeft het einde der bewerking aan enz. 
Men voere de berekeningen ín 5 decimalen uit om een eerste 
benadering der onbekende te bepalen. Het voortzetten der be- 
nadering kan dan door de ín het vorige hoofdstuk geschetste 
impliciete ontwikkeling geschieden. 

7e Voorbeeld. (Zie 4e Voorbeeld). 


Xx— 0.245 3162 sin. x = 27°31’5.23. 


Deze vergelijking is op blz. 50 met behulp der methode van 
Gräffe opgelost. Deze gaf, bij rekening in 5 decimalen als uit- 
komst x, = 35°43/16”. | | 

Als aanvangsbenadering kiezen wij nu x,=835°48’ om de 
berekeningen bij de voorzetting der benadering. zoo eenvoudig 
mogelijk in te richten. 

_Graphisch met behulp van de 3 systemen rechte lijnen vonden 
wij vroeger reeds zeer gemakkelijk en zeer vlug: x = 36°. 

Uit het voorafgegane onderzoek moet nu m.í. de volgende 
conclusie worden getrokken. 

De voordeeligste methode voor de bepaling eener aanvangs- 
waarde van „den” wortel der Kepler'sche vergelijking is de 
graphische methode met behulp van de sinusoïde of de 3 systemen 
rechte lijnen Y) of de numerische benadering door toepassing van 
de formule, die uit de regula falsi is afgeleid. | 


l) Zoolang het cycloïden-tableau niet is verwezenlijkt! 
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Alle andere tot nu toe gebruikelijke methoden komen niet in 
vergelijking met de hier genoemde. 

De methode van Gräffe bewijst hier hare algemeene toepas- 
baarheid, maar vereischt hier toevallig meer becijferingen dan 
de andere. Trouwens bij een vergelijking van eenvoudigen vorm 
zal in het algemeen het in rekening brengen van bijzondere 
eigenaardigheden tot een afwijking van de algemeene oplossing 
leiden. 


C. De voortzetting der benadering. 


Wij stellen: 
fx) =X— esin. x—M=0 
en noemen f= X, —esin. Xx, — M; 
dodmile= € COS. Kon == ESiN. Xp: 
Toegepast zal nu worden formule (21) van hoofdstuk III: 
B f Ni 
X Ie Xo EF IE hi FOREL enne fak ete 
Bij het le voorbeeld: 
Xx — 0.367 5224 sin. x = 27°13/14”.72 

vonden wij vroeger graphisch (met de sinusoïde) x, = 41°. Dan 
wordt f —= — 0.000 6232; f’—= 0.722 6273; f'’=0.241 1164. 

Dus zal: 


0.000 6232 0.241 1164 X 0.000 6232? 
x= 0.715 5850 + G755 6275 9 X 0.722 6273? 


—= 0.715 5850 + 0.000 8624 1 — 0.000 000 1 24 .... 

— 0.716 4473 — 41°2/57.88. | 

Contrôle. 0.716 4473 — e sin. 0.716 4473 — M = 

| 0.716 4473 — 0.716 4473 = 0.000 0000 .... 
Bij het 4e voorbeeld: 
Xx — 0.245 3162 sin. x = 27°31’5”.23 
gaf de graphische methode x,= 36° (3 systemen rechte lijnen) en 
de numerische door achtereenvolgende benaderingen x= 35° 39’. 
De voortzetting der benadering geeft hier bij x, = 35°39’: 
x= 0.622 2098 + 0.001 31164 — 0.000 0001 54....= 
ZD 10 A 400002, 
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Contrôle. Ook 0.623 5213 — e sin. 0.623 5213 — M 
is tot in de laatste decimaal —= 0.000 0000 .... 
Bij ‘het 6e Voorbeeld: 
| X — 0.413 2767 sin. x = 126°47/7.13 
gaf de regula falsi x, == 2.46589. 
WIEZE Pld Le Ie 
De implicite ontwikkeling geeft dan: 


AE 0.005 4538 _ 0.258 4920 X 0005 4538? 
met 0 RET 2 X 1.322 45853 Di 


2.469 9817 = 141°31/10”.29. 


In alle behandelde algebraïsche en transcendente vergelijkingen 
is het groote nut gebleken van de implicite ontwikkelingsformule. 
Bij geen der oplossingen zou Newton's benaderingsformule terstond 
een bruikbare benadering hebben. verschaft; de andere geeft 
echter telkens een resultaat in 7 decimalen nauwkeurig door 
twee correctie-termen. 





V. RECHTSTREEKSCH BEWIJS VAN DE GRONDSTELLING 
VAN DE THEORIE DER ALGEBRAÏSCHE VERGELIJKINGEN. 


le Bewijs. In het navolgende bewijs worden eerst de coëfficienten 
alle reëel gedacht. Het geval van complexe coëfficienten wordt 
daarna tot het eerste herleid. 

De te bewijzen stelling is: 

Elke vergelijking van willekeurigen graad heeft minstens één 
(reëelen of complexen) wortel. 

Wij nemen tot uitgangspunt de 

Hulpstelling: Elke vergelijking van oneven graad (met reëele 
coëfficienten) heeft minstens één reëelen wortel. 

De waarheid dezer hulpstelling blijkt zeer eenvoudig. 

Men kan het eerste lid van zulk een vergelijking altijd positief 
maken door aan de onbekende een genoegzaam groote positieve 
waarde + a te geven; maar het eerste lid kan ook altijd negatief 
gemaakt worden door aan de onbekende een negatieve waarde 
— b toe te kennen, als de volstrekte waarde daarvan slechts 
groot genoeg genomen wordt. 

Het eerste lid moet echter om van een positieve waarde naar 
een negatieve over te gaan (of omgekeerd) zeker door nul gaan; 
de vergelijking zal dus zeker een wortel hebben, gelegen tusschen 
+ a en — hb. 

De algemeene stelling behoeft dus slechts aangetoond te worden 
voor een vergelijking van evenen graad. 

Wij zullen dus trachten aan te toonen, dat aan de vergelijking: 
F(x) = Xx" + Aj trl + Agt? H.H AoriX + Aan= 0... (1) 
kan worden voldaan door een waarde van x, voor te stellen door: 

| rele Ae B (VD 

In de laatste uitdrukkking stellen u en v reëele (positieve of 
negatieve) getallen voor; de uitdrukking is dus in het algemeen 
complex, maar kan voor v =0 toch ook een reëelen wortel 
opleveren. 
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Nu ís: 
F(u tiu) =F(u) + iuF!(u) — Pa) SF. la Flu). (3) 


Aan de vergelijking BE —= 0 kan nu voldaan worden, 
wanneer ies 











u? + (2n) 
Pa) oF 5 zE, Fa) on W= =0. (4) 
en 
OO Oe a 
pn? [2n—1] 
ze 20 


De eerste leden van deze beide vergelijkingen zijn functies van v?; 
wanneer een waarde v = v, voldoet, is er tevens een waarde: 
v=— v,. Hieruit volgt terstond, dat, wanneer een vergelijking 
van evenen graad (met reëele coëfficienten) een oplossing 
X=u+tiv heeft, zij daarnevens ook de oplossing u — iv moet 
bezitten. 

Wij vervangen nu v? door v’ en brengen (4) en (5) in de 
volgende gedaante: 


Ut Yo" EY 0E U Ur 0 (44) 
en | À 
Uv LU Ue Pd Una Ut OE 


Een enkele blik op (4) en (5) zal gemakkelijk doen inzien, 
dat in het algemeen Ux een functie van den Ken graad in u 
voorstelt. 

Tusschen 4a en 5a zullen wij nu v’ elimineeren door middel 
van de interessante methode van Sylvester. 

Hiertoe wordt (4a) achtereenvolgens vermenigvuldigd (gedacht) 
met ov”, o77$ .... tot en met v’. Met (44) zelf meegerekend 
ontstaan aldus n — 2 +1 =n—l vergelijkingen. 

Verder wordt (5a) achtereenvolgens vermenigvuldigd (gedacht) 
metalen Oan tene Samen met (5a) ontstaan 
daardoor dan z Poen SER 

Te zamen zijn er dan 2n — 1 vergelijkingen, die allen volledig, 
afdalend, volgens vi, enz. worden gerangschikt. De 
volledige rangschikking vereischt, dat alle machten van ©’ (van 
ov” tot en met Vv’) in alle vergelijkingen optreden; men be- 
hoeft daartoe slechts nullen als coëfficienten te nemen. 
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De eerste vergelijking van (5a) zou bijv. worden: 

BRU Er Vara OTR: +0 =0. 
De volgende zou dan zijn: 

OP U od + Un Unie!" HOH. H0=0. 
De laatste zou er aldus uitzien: 

ov EO HU 00 HU Un 1=0. 


De vergelijkingen, die uit (4a) ontstaan, zien er geheel 
analoog uit. | 

In het compleete stel van 2n — l vergelijkingen zullen nu 
BD, . tot en met v/ als afzonderlijke onbekenden 
worden beschouwd. Ten opzichte van deze 2n — 2 onbekenden 
zijn de 2n — 1 vergelijkingen lineair. 

Hieruit volgt dan onmiddellijk, dat aan onderstaanden deter- 
minant voldaan zal moeten worden. 





ei OER rds el EO as 0 
0 Kn SP dean 4e Elan 0 0 0 

BOE Ue. 0 0 
et at toeten Wet eltant ea ie maralter se Fadi ddie el est te Wort et) 
hk lj 5: Toten ee Un 0 0 0 Nt 
ONO Ualse0 0 
EE Ea ee U», 0 
eee te We wiel Miete velt u Werle on bed. le, '®' oe vds eve sie ee’ | 
le dr ENE Us, | 


Van deze vergelijking hebben wij slechts den graad noodig. 
Om dezen te bepalen behoeven wij slechts te bedenken, dat (54) 
voor de eliminatie n vergelijkingen heeft moeten leveren en (44) 
slechts n — 1 vergelijkingen. De diagonaal van den determinant 
bestaat dus uit n factoren U, en uit n — 1 factoren U2,; dit 
geeft dus een vorm van den graad: 


nt (n— 1l)2n = In? — n= n(2n — 1). 


Geen enkele andere term van den determinant kan dezen graad 
te boven gaan; een eenvoudige verificatie doet trouwens zien, 
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dat alle termen van denzelfden graad zijn; bijvoorbeeld: 
Ur. U‚”-! heeft den graad n(2n — 1) enz. 5). 

Deze witdrukking voor den graad der vergelijking in u schijnt 
op het eerste gezicht eenigszins ontmoedigend, omdat het onder- 
zoek eener vergelijking van den 2nen graad leidt tot een andere 
van den graad n(2n — 1), die in het algemeen belangrijk hooger 
is dan de eerste. 

Een enkele eenvoudige opmerking geeft hier echter onmiddellijk 
de oplossing van het probleem. 

De factor 2n — 1 van de uitdrukking n(2n — 1) is altijd oneven; 
de factor n heeft altijd een factor 2 minder dan 2n. Hetzelfde is 
het geval met den geheelen vorm n(2n — 1). 

Wanneer wij dus waren uitgegaan van een vergelijking in x, 
welker graad even was, maar slechts één enkelen factor 2 bezat, 
zouden wij ter berekening van u. een vergelijking van oneven 
graad verkregen hebben, zoodat het bestaan van een reëele waarde 
van U zeker aangetoond was. Uit (4a) zou verder het bestaan 
van een reëele waarde voor v? weer volgen. Aan (4a) en (54) 
zou dus gezamenlijk kunnen worden voldaan, dus ook aan 
(4) en (5) en dientengevolge aan de vergelijking (van evenen 
graad) F(x) = 0. 

Wij hebben dus reeds bewezen: Elke evene-machts-vergelijking, 
welker graad slechts één factor 2 bezit, heeft altijd minstens 
één wortel. | | 

Hieruit volgt dan echter onmiddellijk, dat elke vergelijking, 
welker graad 2 factoren 2 bezit ook (minstens) een wortel moet 
hebben. Wanneer men nl. x=u+tiv stelt, zal men ter be- 
rekening van u nog op een vergelijking stuiten, welker graad 
n= n(2n — 1) even is; die graad heeft echter slechts één 
factor 2. Door dus nog eens u = u, + iv, te stellen vindt men 
dan wu, als oplossing van een vergelijking van oneven graad. 

Eveneens wordt v, uit een onevene-machts-vergelijking ge- 
vonden; de berekenbaarheid van wu, dientengevolge die van u 
en ten slotte van x is dus weer aangetoond. 

Op deze wijze voortgaande bewijst men de stelling voor een 
geheel willekeurige hoogere-machts-vergelijking, al moge de 
graad ook een maximum aantal factoren 2 bezitten. 





n—l 


1) Bij gewone, achtereenvolgende eliminatie van vi, v' enz. blijken de coëffi- 


cienten ook steeds gelijkslachtige veeltermen te blijven. 
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Men heeft ín het algemeen te stellen: 
Xi =diurhije: 


Up= Up + iv. 
Door deze uitdrukkingen op te tellen vindt men dan: 
B EE Warte doeh nt IN. 

De waarde van up is reëel; boven merkten wij op, dat daarbij 
dan een reëele waarde van V? behoort. V zal weer òf reëel of 
zuiver-imaginair moeten zijn. Derhalve is bewezen: 

Elke vergelijking van evenen graad heeft altijd een wortel, die 
òf reëel òf complex ís. 

Bij het voorafgaande betoog zijn de coëfficienten der be- 
schouwde vergelijking steeds reëel ondersteld. Bij aanwezigheid 
van complexe coëfficienten behoeft men slechts te stellen: x=p gi; 
om na substitutie het reëele en het imaginaire gedeelte van het 
eerste lid der vergelijking afzonderlijk = O te stellen. Door 
eliminatie van p of q ontstaat dan weer een vergelijking met 
reëele coëfficienten, zoodat men op het vorige geval terug- 

gekomen is. 

2de Bewijs. Men kan zich afvragen of aan de vergelijking: 

OE ANA re He aderen Ant Aar =O 
ook kan worden voldaan door de wortels der vierkants-vergelijking: 
| Xe — PX —gq=0. 

Wij kunnen daartoe stellen: 

BOE et Aj Oep Apster. nt Anikrt Aan = 
Eg ES Pa Ee et Aansartt dan--2): 
Het is terstond duidelijk, dat dit mogelijk is, onder voorwaarde, 
dat aan onderstaande vergelijkingen kan worden voldaan: 
| =p A 
A= Apt A tg 
As = Aap + A; + ag 
of wel jas = asp + À, + 434 
CRN ee 


A —Pp=À 
ap g=A, 
A3 — Ap — Aq = A; 


Qan—2 — An—3P — Aan-4 = Aon 
ardin — Aon] == Aon 
— An = Aon 


Qan-2 = Aan—3P + Àon-2 + Aon-44 
Aan an AU2n—3q Aars 0 
Ano Ein Aon TE 0. 
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Het blijkt onmiddellijk, dat alle onbekenden a,, a, enz. achter- 
eenvolgens kunnen worden uitgedrukt ín p en q en in de gegeven 
coëfficiënten A, A, ....Aor, zoodat men ten slotte twee ver- 
schillende vergelijkingen in p en q verkrijgt. Als het nu mogelijk 
blijkt aan deze beide laatste vergelijkingen te voldoen, zal het 
geheele systeem (A) oplosbaar zijn. Dan is de aanwezigheid van 
den factor x°—px—gqgq en derhalve het bestaan van twee 
wortels van de beschouwde vergelijking aangetoond. 

De graad der opeenvolgende uitdrukkingen van a,, 4, enz. in 
p en q wordt opgegeven in onderstaand tafeltje: 


graad vans. dy. „dw 10 0e ne Mans SR 
in p Faith Aers BiA A 
in q O1 1 2e Nn En nk 


De beide laatste vergelijkingen van het systeem (À) zullen dus 
ten opzichte van g respectievelijk van den „ — len en den nen 


graad zijn. 
Zij kunnen worden voorgesteld door: 
Pgt UF Pegel dst Pand Tian =O 
ens Pag Pag ot eet Ponad Pane en (7) 


(Px is een polynoom in p van den ken graad.) 

De eliminatie van q geschiedt volgens de methode van Sylvester 
op precies dezelfde wijze als die van v’ in het /e bewijs. 

De resulteerende vergelijking voor p wordt: 








PUPS Pasir 50 Li 05 ANO 
OND OP ER Pine FOM 0 
OO ARSEETEED o D Pant 0 0 
RENT Aen 0 De nt piel 0 
BARDO Dn. Pin "OO 0 a 
OM PRS REAR HA: EER Pis a0 0 
OR Oe ARAS Po Pt NE NES ABE Pan 0 
0 Orr ORO de Hie AE LENEN 


De graad van deze vergelijking zal weer = n(2n — 1) blijken 
te zijn, nl. gelijk aan dien van een willekeurigen term er van. 
Omdat (6) in het geheel n vergelijkingen voor de eliminatie heeft 


89 


moeten leveren en (7) in het geheel n — 1 vergelijkingen, zal de 
diagonaal = P, P!' zijn en dus van den graad 


n_ 2n 
nt 2n(n — 1) = n(2n — 1) 
zijn. Elk der overige termen zal van juist denzelfden graad ‘zijn !). 

Men ziet weer terstond, dat de graad van de resulteerende 
vergelijking een factor 2 minder heeft dan die der oorspronkelijke. 
Als de graad der oorspronkelijke vergelijking slechts één factor 
twee bezit, is de resultante in p van oneven graad; er bestaat 
dan een reëele waarde voor p. Uit (7) volgt dan eveneens een 
reëele waarde voor g. De reëele factor Xx? — px — q van het 
eerste lid van F(x)=0 is dan aanwezig en het bestaan van 
twee wortels van deze vergelijking van evenen graad is dan 
bewezen. Deze twee wortels kunnen reëel of complex zijn. 

Wanneer de graad van F(x) = 0 twee factoren 2 bezit, zal de 
resultante in p nog een graad met één factor 2 hebben. Deze 
resultante heeft dan echter een deeler van den tweeden graad; 
er zijn dan twee waarden voor p. Door één er van te substit 
tueeren in (6) en (7) verkrijgt men twee vergelijkingen voor g, 
die op bekende wijze zoover men verkiest, — bijv. tot een ver- 
gelijking van den eersten graad, — te herleiden zijn. Er bestaat 
dus zeker weer een waarde voor q en GE weer een tweede- 
machts-deeler x? — px — g. 

Op boven beschreven wijze kan men de stelling voor steeds 
grooter aantal factoren 2 in den graad aantoonen, zoodat zij 
geheel algemeen bewezen is. 

Opmerking. Het valt zeer gemakkelijk te bewijzen, dat een 
vergelijking (met reëele coëfficienten), die een oplossing x= u + iv 
heeft, ook een wortel x= u — iv moet bezitten ®). Zij zal dan 
in het eerste lid den factor 
| (Xu — iv) — u + iv) = x* — ux + u? + v?, 
dus een reëelen 2e machts-factor moeten bezitten. 


Verificatie en Toepassing der gevonden resultaten. 


Nemen wij tot voorbeeld de 4e machts-vergelijking: 
Ere tende Als sel). 


1) Zie de noot op blz. 86. 
2) Zie het Je bewijs. 
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1°. Wij stellen: Kd (U, 

Volgens de formules (4) en (5) zijn u en v gegeven door: 
us Au? d Azu? + Agu + A, — v(6u? + 3A,u + Aj) duf =0 
en 
Aus + 3A,u? + 2Asu + Ag — U? (4u + Aj) = 0. 

De eliminatie van v° leidt tot een vergelijking van den 6en graad: 
(4u HA) (ut HA, us HAou? + Agu tAa) + (4u? +3A, u? + 2Aou A3)? — 
(4u? + 3A u? HAgu + As) (6u? +3A u + As) (Au HA) =0. 


In het geval, dat A, = 0 ad En de laatste vergelijking over in: 


uê ien gu 1e za u? ee m0 





16 


Deze is juist de hulpvergelijking van de Euler'sche methode. 
2%. Wij kunnen het eerste lid gelijk stellen aan: 


(X° — PX — 0) (XT ajk T 43). 
Volgens de formules van het systeem (À) ontstaan dan: 
(2p + Aj)g + p° + Ap? + Asp + Ag =O 
en 
q°+(3p°+2Aip + As)g + pt + Apt Arp? + Asp tr A, =0. 
De eliminatie van q voert tot de 6e machts-vergelijking: 
p Ay) (pt Aip + App? + Asp Aa) + (p°tArprtAaptAg)— 
— (PP +Aip' + Arp + Aj) (3p* + 2Ap + As) (2p + Aj) =0. 
De graad n(2n — 1) komt uit, want 2(4 — 1) = 6. 
Wanneer A, = 0 is, gaat de vergelijking over in: 
p° + 2A,p* + (As? — 4A)p? — A? =0 
d.w.z. juist in de resolwente van Descartes. 


Afleiding der algemeene eindvergelijkingen bij het lle bewijs. 


Wij zullen ons bedienen van de notatie: 
Falp)= pr Ap" VEA pen An 1D A 
dus Fra1i(p) =p"* UF Ar pir Aap olst Aap Tt An 


waaruit volgt: 


PAD) DE DIAR: e 1E ANT ON (10) 
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Uit (10) kan men afleiden: 


ed Bal p) PE nl p)eheFalp); 0. (II) 
eveneens: BOEDA Ent MO, Mate LZ) 
Uit (12) kan weer afgeleid worden: 
EERDE DIP ni Pin A nor 1e (13) 
in het algemeen: 
B EE (pm (py ets NE 10 (14) 


Wij kunnen nu ee dat een willekeurige onbekende a, 
van systeem (A) aldus kan worden voorgesteld: 


= F‚(p) + qFn—i(p) + Sin xp) + A A REE. 


Fe (P) EEE) 


Deze formule heeft betrekking op het geval, dat k even is; 


wanneer A oneven is, zal de hoogste macht van g gelijk zijn 
En, 
2 
Om (15) te bewijzen zullen wij aantoonen, dat zij waar is voor 
An+1, wanneer zij geldt voor aj: en voor ap. 


Bij (15) voegen wij dus: 
Ani = Fri(p) + qF’n_o(p) + Ce 3(p) AT d ee 4(P) + 


4 el 


eee ED (16) 
Fn Mt 
5 2 





Nu ís volgens systeem (A): 
Gri Dühr Une re Apr: 
Dus wordt a„‚1 gevonden door beide volgende uitdrukkingen 
samen te tellen: 


Aert PFa(p) + paF'n(p) + pi PP) + pr SF pt 











4 
2 
gE iP) HPP rl) + SED Spr). 
tn 2 
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De opvolgende termen der resulteerende formule zullen zijn: 








Ansi t pFn(p) = Frsilp); . …. . (Zie 10) 
bk ‚ … (zie 12) 
PPP) + ral) }= PFP: … Gie 13) 
h h 
2 id KC pi h 
Eler ot dek lo aan zie (14) 
bee 5 2 ij JN) Di. 


Hieruit volgt: 
8 
Ansi = Pre) Ha) Hip NT eee einde 


d.w.z. juist hetgeen te bewijzen was. 

Een eenvoudige verificatie bewijst terstond, dat de formule 
waar is voor h=l en h= 2; daaruit volgt dan de waarheid 
voor alle geheele waarden van 4. 

Wij hebben dus bijvoorbeeld: 


amn-a=Fan-alp)+GF ana) +5 4e zalle en 





Arn-= Fan dp) Haan) + 5 7 Fl) Jp). 


Ee — Lj 
Uit: Aan-2P + Aon-3q + Aon =O (Of aan = 0) 
volgt: 


Parr(@)H0P lp) AF (P) + 








n—l 
RI Bon Er "wo EN hbe 
Verder kan Aan-2q + Aon = 0 


vervangen worden door: 


QF an 0) + LF 3(P) + LOA 


trent 0 + A= 0 A 
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Eindelijk geeft p X (17) + (18) ons de vergelijking: 
ie 
Fan (p) + Paar (P) + jg Fia) + 





rd — Frl(p)=0, . . (19) 


die men onmiddellijk vinden kan door te bedenken, dat a1, = 0. 
De eindvergelijkingen zijn volgens hetgeen voorafgaat (17) en (19). 
Zij geven rechtstreeks de bijzondere vergelijkingen, waardoor 
wij de hulpvergelijkingen van Descartes en Euler hebben afgeleid. 


VI. NIEUWE RECHTSTREEKSCHE OPLOSSING VAN DE 
CUBISCHE VERGELIJKING. 


De bekende oplossing van (CARDANUS gaat uit van den 
bijzonderen vorm: 


NTP Oe. 
De onderstaande methode neemt als uitgangspunt de vergelijking: 
ttid A NE oh BBA AD 0 … Aten 


De algemeene derde-machts-vergelijking is door een lineaire 
substitutie altijd tot bovenstaande gedaante te herleiden. 

Zij de algemeene vergelijking bijvoorbeeld 3 + ay? + by + c=0. 
Door de substitutie y — Xx + p- ontstaat: 

Gesld)Betal kp) DOREN 

of X8 + (3p da)? + (Bp? apa tp tapt bp 

Om den vorm (1) te krijgen moet nu: 

(Sp 0){3p2 H-2ap FOCON 

Bij uitwerking vallen de termen met p? en p? uit. Er blijft 

alleen over: 
2(a° — 3b)p + ab — Ic = 0, 


waaruit volgt: 


OL AAO 
HA Sat 285) A 
In (1) stellen wij nu: | 
en DS VB (+ of — naar kêuze). re 


Er ontstaat dan (na verdrijving der breuken) een cubische 
vergelijking in y, waarin het drievoud van den coëfficient van y 
gelijk is aan het quadraat van dien van y?, zoodat het eerste 
lid in den vorm (y + a)? te brengen is. 
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De afleiding van die vergelijking is zeer eenvoudig: 
slot vB} ot + P +878 
EV, dr jy y 
6A B 


3Ax® BAB 
y y 


9Bx = 7 + 3BV/B 
AB = AB 


0 RL SGEE 6V B(A+/ B) 
y y 





EE + 4B(A + /B) 





3 3 1 
5) 2 De rm 
a VBV Tat Ba rvB 
Hieruit volgt: 
| il je EE I 
JT 9yBl8BVB 4B(A + VB) 
ee Pr ER 
ki Bn er U, 
Deze betrekking geeft terstond drie waarden voor y: 
1 1e AVB 1 | Bed 
0 pe ì NE ebde RENS 
U A toyBrovs VaryB Yaya tt Varvs) 
AVENE B 1 
0 EE Pe ë 9 vnd al AN EEL 
LTE ov VAryB INT 575 
vd tse dE 
AVB A 4D! 

















(4) 


l 1 +i/3 A—VB 
0 Egli, hk, 
OVB 4yB VATVE 
Uit (2) volgen dan de 3 waarden voor x. 
Voorbeeld: 





Ket Dl 93 103 = 0; 
Del j Pp en neem Le A dk EI 2 
Td PS T4A1F3X08) T 1440 
xH(3p 2 + (Bp? + 42p — 3)! + p3 + 2Up? — 3p +103 =0 
gaat dan over in: | 
| xt TK +3 +9 =0. 
Nid en Bl. 
Uit (4) volgt dan: 


DL 0 Tof EDI 05 
Yr ao or os " Ar " ‚JJ. 





Ì 


96 


Met behulp van (2) heeft men dan: 


LE 














Rr VOS 05 TT l, zoodat #, = MIE OË 
Dus zal: 
pe SVOD, 0 5 VAO DEN 6) 
LVO O5 ODE DAN 
Verder is: 
Lif 373 BAT EE (0 NN 
jee Se Te de d Le : 
OR 1 10 DVO 
N= Hl— i/3) VO = EN 


zoodat: %= 3 — Ge Oje EL 
k _1F(l—i3VOIL 1H(l if) VO1 

Op dezelfde wijze blijkt: 

Bs. 1 + 3(1 + i/3) 0.1 D 
Br lt (lot VS OL enen Es 

Er is zeer veel voor om de algemeene 3e machts- VE 
in den bijzonderen vorm (1) te brengen. 

Vooreerst geldt ten opzichte van de bestaanbaarheid of onbe- 
staanbaarheid der wortels de volgende, zeer eenvoudige, 
Stelling. Is B positief, dan zijn twee worteis complex; 

is B negatief, dan zijn alle drie wortels reëel; 
is B —= 0, dan zijn ook twee wortels = 0 en de derde reëel. 

Bewijs. Het eerste deel dezer stelling volgt terstond uit de 

beschouwing van (4) en (2). De oplossing y, zal, als B positief 


(6) 








is, reëel zijn en dientengevolge ook Xx” = 5 + VB en ten stottesn 
1 


Zoowel y, als ys en dientengevolge weer X en X, zullen echter 


complex zijn. 
Het tweede deel der stelling blijkt door x zchtecc en 


gelijk te stellen aan O, — 3A, + V—3B en — V— 3B. 


Dit geeft: 
MOOT Oek wordt de waarde van het le lid van (1) = + AB 
„ X=—3Â vs om end ROR eenn 
„nam 3Bj gd onde ee en a 
Xb HEAD eten A an A 


n 


Er moet dus een reëele wortel liggen tusschen O en + V —3B en 
ook tusschen O en — V—3B. Er zijn dus 2, maar dan ook 3 reëele 
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wortels. Van de beide eerste is de eene positief, de andere 
negatief, het teeken van den laatsten hangt af van dat van A. 

Wij kunnen nu aan bovengenoemde stelling nog toevoegen de 
Stelling: Is B negatief, dan heeft (l) een reëelen wortel 

a) tusschen 0 en — 8A; 

b) tusschen O en + V— 3B; 

c) tusschen O en — V— 3B. | 

De gevallen a en b evenals de gevallen a en c kunnen den- 
zelfden reëelen wortel aanduiden, maar de aanwezigheid van 
3 reëele wortels volgt vast uit b en c. 

Toepassing: KEO d leet kli): 

DIE Heet p en teem p= Giet 658 

Verminder dus de wortels der gegeven vergelijking met 5.5 en 
hiertoe eerst met 5. Uit de coëfficientenrij 

1 6 11 1 

is terstond duidelijk, dat het teeken van 3B altijd positief zal 
blijven: de vergelijking heeft dus 2 onbestaanbare wortels. In 
het algemeen behoeft men slechts het aantal geheelen van p van 
de wortels af te trekken om terstond te kunnen zien, wat het 
teeken van 3B zal zijn. 


Nieuwe Goniometrische Oplossing. 


Aan de gedaante (1) sluit zich verder een zeer eenvoudige 
goniometrische oplossing aan: 
Beschouwt men de formule: 


tg. 3x — Stg. 3xtg. °x — 3tg.x + tg.3x=0 ... (8) 
in verband met de algemeene cubische vergelijking: 
EEN an 1e Ona le EERE (9) 


dan blijkt (8) identiek te zijn aan (9) op voorwaarde, dat men 
mag stellen: 

a=——3tg. 3x, b=—38 en cC=tg.3 >. +10) 

Het gelijktijdig bestaan der 3 betrekkingen (10) vereischt, dat 
tegelijkertijd: 


LE OE OTU DT dte el dien (11) 
Aan (11) wordt ook voldaan, wanneer tegelijkertijd 
Ag BL Din rehans an te de rg nt sdan e (12) 


2 
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De eerste dezer beide betrekkingen is juist vervuld in onzen 
vorm (Ì); de tweede kan zeer gemakkelijk daaraan toegevoegd 
worden. / 

De goniometrische oplossing wordt dan deze: 

Breng de algemeene cubische vergelijking in den vorm: 


KS HTB ANE TGD HE ARTE 
Stel dan x= py, zoodat 
JAN 3B AB 
dn fn on eer al 
d p y p” p° 
Stel Ee =—d en dus p= tV—B (+ of — naar keuze). 
P IN 
De vergelijking in y wordt dan: 
‚JA Eee 
3 LEREN KT PE Dt erk 0, 
JASON EN 0. 


Stelt men nu y =tg.y’, dan volgt uit (8), dat: 


ARL RE) A abs —— Ee o 

ter av GLR: zoodat y — g Pgtg. F VB rn 
en ten slotte: 
Emil Íb tg. + Ed d tx 60° (13) 
y=ten belg. F rp + Bd | 
Men heeft voor x dan nog x= py, dus 
nh 1 AAR k 
X= + VB te jg bete. + VR + A X 60 , 


Uit (13) volgt terstond: 

1°. dat deze oplossing geldt voor het geval, dat B negatief 
is, dus dat de 3 wortels reëel zijn en derhalve voor het z.g. 
onherleidbare geval; 

20, dat er 3, maar ook niet meer dan 3 reëele oplossingen 


gevonden worden. Als men het teeken in + ve zoo gekozen 


heeft, dat deze vorm en dus ook 7 bete: + B positief is, 


| Krt HA | 
heeft men naast tg. E bgtg. V=B nog de waarden: 
® 1 rr Ä | l _—_—— JN le) 


Alle andere oplossingen zijn tot een van deze drie te herleiden. 
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Voorbeeld: 
X8 — 30x? + 96x — 48 = 0. 
Stel Pasen ot 0 
zoodat 
x°H(3p — 30)!” + (Bp? — 60p + 96)X’ + p3 — 30p? + 6p — 48 =O 
en doe GD Pa Aero dOr eg 
(a? — 36) 2(900 — 288) 
xj — Mr — 1 +32 =0. 
Nu is dus A=—8 en B=—4. 


De oplossing is: 
X =2+2tg.$tbgtg. 44; 
Xa = 2 + 21e. Ht bgtg. 4 + 60°4; 
X= 2 + 2tg.$$bgtg. 4 + 120° 4 
Dit geeft: 
Xx = 2.9463035; x, — 2643746; x, — 0.616232. 


Contrôle. De som der gevonden wortel is 29.99999 (55); 
deze ligt zeer nabij de gegeven waarde: 30.00000 00. 





VII. ONDERZOEK VAN DE RECHTSTREEKSCHE OPLOS- 
BAARHEID VAN HOOGERE-MACHTS-VERGELIJKINGEN 
DOOR MIDDEL VAN DETERMINANTEN. 


Stelling. De eliminatie van z tusschen G(z) = 0 en H(z) = y 
levert altijd een resultante R(y)=0, die van denzelfden graad 
is als G(z) = 0. 

G(z) en H(z) stellen geheele functies van z voor, waarvan 
alleen G(z) van eindigen graad behoeft te zijn. 

Het bewijs is eenvoudig. De nieuwe onbekende is een één- 
waardige functie van z; voor elk der wortels van G(z) = 0 kan y 
dus slechts een enkele waarde hebben. R(y) = 0 kan dus slechts 
n wortels hebben, als G(z)=0 die heeft en moet dus met 
G(z2) = 0 van denzelfden graad zijn. 

Opmerking. Men kan de stelling uitbreiden tot: 


De eliminatie van z tusschen G(z) = 0 en y = Ne moet een 
resultante opleveren R(y) = 0, die van denzelfden graad zal zijn 
als G(z) = 0. 

Toepassing. Combineert men: 

zit Ast lede ta A  et AeR 
met ela; e 2 ese he Dao al EB 


(of in het algemeen met z"-? 4 Bj ar Pel th ont Biopic 
dan zal R(y)= 0 van den nen graad zijn. 

De coëfficienten van-R(y) = O kan men beschouwen als functies 
van de, voorloopig onbepaald gelaten, n — 1 grootheden B,, 
Deen. 

Men kan nu trachten over B, Bs, .... Br zoodanig te 
beschikken, dat de coëfficienten van yr-l, pr? ....yalle=0 
worden. Hiertoe heeft men n — 1 vergelijkingen met evenveel 
onbekenden op te lossen. Daarna zou men dan y vinden uit een 
binomische vergelijking van den nen graad en ten slotte ter be- 
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rekening van z een vergelijking van den n — len graad over- 
houden. De oplossing eener vergelijking van den ren graad zou 
op deze wijze teruggebracht worden tot die eener vergelijking 
van den „— len graad. 

Werkelijk gelukt deze methode bij vergelijkingen van den 
tweeden, derden en vierden graad. 


L. De vierkantsvergelijking. 


De resultante der eliminatie van z tusschen 
sal —= 0 en 
dBm. 0! IS 
ie RE icrÂ, Às En 
1 (B, —y) O0 O0 of (B, —y)’ —A,(B, —y) + A, = 0. 
| 0 bonden) | 
Uitgewerkt, wordt zij: 
yv? —(2B, — Aj)y + B,° — AB, + A, = 0. 
Men heeft nu te stellen: 


B, — A, =0 of B, =S 


Voor y vindt men de binomische vergelijking: 
2 Eide 
pe En — As, dus y= + 4VA,? — 4A, 
Uit z + B, — y = 0 volgt dan: 
A 


de Ae UA, 


1. De derde-machts-vergelijking. 


Bij ZAAR AR oi Ans 0 
en zere 0D 
behoort als resultante: 

1 A, Às As 


Ged mad eh hed 
mw 
Te 
es, 
| 
tE, 

(en; 
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Het is met een oogopslag duidelijk, dat bij uitwerking y° een 
coëfficient zal hebben, die in B, en B, Lineair is, terwijl de 
coëfficient van y een guadratische functie zal zijn van B, en Bs. 
Men kan dus door oplossing eener wvierkantsvergelijking de termen 


met py? en y doen verdwijnen en dus y oplossen uit een binomische 
vergelijking y? = C. | 


IL. De vierde-machts-vergelijking. 


De resultante der eliminatie van z tusschen 
ze JAE TA tT ARe Tt An 
en 
23 J-Bie nt Bien Dit Vid 
is de onderstaande determinant: 


BONE PO MP RLNDOE As 0 op 

0 1 ARE bi A, A, ol 

0 0 BAG A, ACME 

LGB: BAB 0e 0 0 0 RED 
0 1 Blin HiB BE Ren 0 

0 0 1 B, Be (BES R0 

0 0 0 1 B, B, (B;-—) 





Wilde men hier nu de coëfficienten van y°, y? en y = 0 maken, 
dan zou men voor B,, B, en B; drie vergelijkingen hebben op 
te lossen, die respectievelijk van den len, Zen en 3en graad zijn. 
Zulk een stelsel zou ín het algemeen tot de oplossing eener 
vergelijking van den 6en graad voeren. 

Men kan de methode echter nog blijven toepassen, als men 
slechts verlangt, dat in R(y) = 0 de termen met py en y zullen 
ontbreken, zoodat R(y) = 0 en aldus zal uitzien: 


VAG RON = () 
en nog als een vierkantsvergelijking kan worden opgelost. 


Er behoeven dan ook slechts 2 onbekende coëfficienten te 
worden ingevoerd, waardoor nog een vereenvoudiging bereikt wordt. 


Bij Zhi Agee ele eenn Anr (0) 
neme men dan: 
Bie + (Bet) 0; 


De resultante wordt: 


RR A, A, A, 0 

RES A, A, A UNA. 

l Ber ant De ss) Ò Ò 

1 BB 0 0 mast 
RE 40 1 Bee 

0 0 LRAD (Be) 


De coëfficient van y* zal weer in B, en B, lineair zijn, die 
van y is van die grootheden een functie van den Sen graad. 
B, en B, zijn dus te berekenen en derhalve is de vergelijking 
van den 4en graad rechtstreeks op te lossen. 

Bij de vijfde-machts-vergelijking, — en ook bij vergelijkingen 
van hoogeren graad, — blijft de methode in gebreke. Geen 
wonder: reeds door AgerL is bewezen, dat de rechtstreeksche 
oplossing van vergelijkingen van den vijfden en hoogeren graad 
in het algemeen onmogelijk is. 

Het komt mij voor, dat dit feit met meer vertrouwen door een 
„beoefenaar van de hoogere algebra zal worden aanvaard na lezing 
van boven gegeven algemeene beschouwing over oplosbaarheid. 
Abel's bewijs is verre van elementair en in alle, mij bekende, 
leerboeken wordt tot nu toe zijn stelling steeds plotseling zonder 
eenig bewijs gememoreerd nadat de zeer bijzondere methoden 
van Cardanus, Euler en Descartes zijn vertoond. 
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Aanhangsel A. 
GRAFIEK DER SINUSOIDE 


ter (benaderde) graphische oplossing van de Kepler'sche Vergelijking. 
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Aanhangsel B. 
GEWIJZIGDE SINUSOIDE 


ter (benaderde) graphische oplossing van de Kepler'sche Vergelijking. 
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Aanhangsel C. 
DRIE SYSTEMEN RECHTE LIJNEN 


ter (benaderde) graphische oplossing van de Kepler'’sche Vergelijking. 











dd: HH ze 








/ Es, ne Lm: 
EDE EL E 4 


(RES EE AE Ad AE HE 











mm’, 


HH OD 


En 
BERNER EUGEDEE 
FHT dE 

EEDE 


/ 
Zi 
mm: 
E 
LJ 
LÂ 
4 
4e 
LL 


EEE: 


Pas cinrdrHaas 


EHH 
EH 
ps H 


EED SA 


EERDERE SE 


on B 


mans 
ann 


EH Ö m 
/L FHA 
Ä EEH 


FH 
ann rinv.ad, / 


Ps 
TST 


Ehle) 


HA ian mens " 
J mars ware ue 
e Z| Zi m 
LH | B7,m m 
auv/nv ann Lj T 
Darm AB/INDE , 
/ DORY, wrm ms 


de 


/ CH W. 
mi 








7 


vir 


,# 
k 
8 

E 


wipe vree 








ed warner dik in 


rvs wopr nn en ml 1 s. 


er ee ed 








